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O : 1 INTRODUCTION 

^ . The heavy electron state in intermetallic lanthanide and actinide compounds has its origin in 
^ '. the interplay of strong Coulomb repulsion in the 4f- and 5f- shells and their hybridisation with 
; conduction band states. Heavy quasiparticles have mostly been observed in Ce- and U-based 
Q I intermetallic s which are therefore at the focus of this review. Ideally the ground state in these 
O ■ strongly correlated electron compounds may be described as a Landau Fermi liquid (LFL) state 
>• with large enhancement of the effective mass and associated large values of the linear spe- 
^ cific heat coeffient 7, Pauli susceptibility Xo and A-coefficient of the resistivity. It has recently 
^ ! become clear that the microscopic origin of mass enhancement is quite different in Ce- and 
" " ' U-compounds, as described by the Kondo lattice model and dual model respectively. 

The renormalised LFL state is however prone to instabilities at low temperatures due to 
residual heavy quasiparticle interactions and the sharpness of the Fermi distribution. The most 
common instabilities are the spontaneous appearance of (spin-)density wave (SDW) and super- 
conducting (SC) order parameters which break at least some of the underlying spatial or internal 
symmetries like time reversal or gauge symmetry. Via their associated gap functions both type 
of order parameters are reflected in the modified energies and densities of low lying excitations. 
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This also changes low temperature properties in the ordered state and may in turn be used as 
a means to obtain information on the symmetry class of the order parameter. Both SC and 
SDW gap functions A(k) may be conventional or unconventional depending on whether their 
k-dependence has the same or a lower spatial symmetry as the Fermi surface. In the latter case 
particles (SC) or particles and holes (SDW) pair preferentially at neighboring sites and avoid 
the strong on-site repulsion of quasiparticles. For this reason the SC states in Ce- and U-based 
HF compounds are usually of the unconventional type. Their intriguing properties have lead 
to a rich and flourishing field of reserarch, not only in the genuine HF compounds but also in 
oxides, ruthenates and organic solids. 

The stability of the ordered phase may be influenced by changing microscopic control pa- 
rameters via application of pressure or by chemical substitution. In this way the SDW state may 
be tuned to a magnetic quantum critical point (QCP) where the staggered moment vanishes. 
On theoretical grounds it has long been suspected that unconventional SC pairing is favored by 
magnetic instabilities, triplet pairing in the ferromagnetic and singlet pairing in the antiferro- 
magnetic case. This correlation, however, cannot always be upheld for real HF materials. 

The accumulated evidence in the Cel22 (e.g. CeCu2Si2), Cell5 (e.g. CeCoIns) and Ce218 
(e.g. Ce2CoIn8) classes of superconductors seems to vindicate this picture, since the SC appears 
mostly in small 'domes' around the QCP. The connection between critical spin fluctuation prop- 
erties and the stability of SC pair states has been analysed within strong coupling Eliashberg 
type theories. The QCP has additional implications beyond SC. In the normal state above the 
SC dome thermodynamic and transport behaviour show distinctly anomalous non-Fermi liquid 
(NFL) behaviour as function of temperature and field. This is thought to be the result of a dress- 
ing of quasiparticles with soft spin fluctuations leading to entirely different scaling exponents 
for specific heat, resistivity etc. as compared to the LFL state. This scenario is accepted for 
most NFL anomalies in Ce-HF compounds. Alternatively they may be caused by the existence 
of a 'pseudo-gap' associated with a 'hidden order' parameter, e.g. an unconventional SDW. This 
may indeed play a role in CeCoIns and Celrlns in a similar way as invoked for the pseudo-gap 
phase of underdoped cuprates. 

In U-based HF superconductors quantum critical behaviour does not play an important role 
with the possible exception of UBei3. Instead of appearing close to the destruction of SDW 
order, SC in U-HF compounds is mostly embedded within a stable AF phase of reduced (some- 
times very small) ordered moments. Like the different origin of mass enhancement this may be 
connected with the multi-orbital structure and multiple occupation of 5f shells. Instead of the 
Kondo lattice picture, a dual model with partly localised and partly itinerant electrons caused by 
a strongly orbital dependent effective hybridisation is more appropriate for U-HF compounds. 
In turn this suggests a new SC mechanism: Pair formation is caused by the exchange of mag- 
netic excitons which are CEF excitations of the localised 5f-electrons that have acquired dis- 
persion due to intersite-exchange. Contrary to the critical spin fluctuations in Ce-compounds 
these propagating bosonic modes are not overdamped. This mechanism has been vindicated in 
UPd2Al3 by the complementary results of quasiparticle tunneling and inelastic neutron scatter- 
ing (INS) experiments. An equivalent experimental support for the spin fluctuation mechanism 
is still lacking. 

This article reviews the present understanding of Ce-based and some of the U-based HF 
superconductors, with emphasis on the former. We will not discuss the physics of HF su- 
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perconductors containing other 4f or 5f elements like Pr skutterudites or trans-uranium based 
superconductors. For Ce-compounds we focus on the connection to magnetic quantum critical 
behaviour and the implications of the dual 5f-electron nature are discussed extensively for U- 
compounds. In sect. 2 we give the theoretical foundation of the different heavy quasiparticle 
origin in Ce- and U- compounds and present a brief summary of NFL properties in the normal 
state close to a QCP. In sect. 3 we discuss the known classes of Ce-HF superconductors with 
emphasis on the Cel22 compounds. We also discuss the new non-centrosymmetric SC CePtsSi. 
In sect. 4 we give a description of the different microscopic pairing mechanisms present in the 
Ce- and U- intermetallics. In sect. 5 we discuss only those two U-based HF superconductors 
where new results have been obtained recently, namely UBei3 and UPd2Al3. Finally sect. 6 
gives the summary. 

Of the many already existing review articles and monographies on HF systems and their 
superconducting state we mention here a small selection: The normal state properties of HF 
metals are at the focus in [1-4]. General reviews on HF superconductors are given in [5, 6] 
and the monography [7]. The exceptional case of the multicomponent SC UPt3 is discussed in 
detail in [8, 9]. 

2 THEORffiS FOR THE NORMAL HF STATE 

The theoretical understanding of superconductivity and magnetism in the heavy fermion sys- 
tems is still in the state of rather schematic or illustrative models. A major difficulty is that the 
normal state quasiparticles can sofar be described only within effective single-particle renor- 
malised band pictures with empirical input parameters. For some compounds like UBei3 and 
Ce- compounds close to the quantum critical point the SC transition may even take place in a 
state where the low energy excitations are not simple LFL quasiparticles but are dressed by soft 
spin fluctuations. This is witnessed by the observation of non-Fermi liquid (NFL) behaviour 
in thermodynamic and transport quantities. A fully microscopic description in the case where 
inter-site effects become important is not available. Typical lattice effects are the formation of 
coherent heavy quasiparticle bands whose Fermi surfaces were observed experimentally. We 
first discuss the renormalised band theory which provides a way to describe the coherent heavy 
quasiparticle bands within a Fermi liquid approach. The latter can be calculated from single- 
particle Hamiltonians where the effective potential is constructed to account for many-body 
effects. The residual interaction among the quasiparticles eventually leads to the instability of 
the normal Fermi liquid phase. 

2. 1 RENORMALISED BAND THEORY FOR HF Ce-COMPOUNDS 

The Landau theory assumes that there exists a one-to-one correspondence between the states of 
the complex interacting system and those of a gas of independent fermions, which may move 
in an external potential [10-12], [13]. The single-particle orbitals and energies are determined 
from an effective Hamiltonian. The characteristic properties of a system are reflected in an 
effective and not necessarily local potential which describes the field of the nuclei and the 
modifications arising from the presence of the other electrons. The essential many-body aspects 



3 



of the problem are then contained in the prescription for constructing the effective potentials 
which have to be determined for specific problems. 

The Landau theory of Fermi liquids is a phenomenological theory. The characteristic prop- 
erties of the quasiparticles which can hardly be calculated from microscopic theories are ex- 
pressed in terms of parameters which are determined from experiment. The quasiparticle ener- 
gies in a crystal 

e(k) = v;'(k)-(k-k^(k)) (1) 

are given in terms of the (anisotropic) Fermi wave vector kp and the Fermi velocity yp which 
depend upon the direction k. The key idea of the renormaled band method is to determine these 
quantities by computing the band structure for a given effective potential which accounts for 
the many-body effects. The periodic potential leads to multiple- scattering processes involv- 
ing scattering off the individual centers as well as the propagation between the centers. The 
characteristic properties of a given material enter through the information about single center 
scattering. They can be expressed in terms of properly chosen set of phase shifts {T]y (e)} spec- 
ifying the change in phase of a wave incident on site i with energy E and symmetry v with 
respect to the scattering center. Within the scattering formulation of the band structure prob- 
lem the values of the phase shifts at the Fermi energy {r}'y{eF)} together with their derivatives 

I [drfy I de) | determine the Fermi wave vectors kf and the Fermi velocity \f- The scattering 
formulation of the effective band structure problem provides a highly efficient parametrisation 
scheme for the quasiparticles. To further reduce the number of phenomenological parameters 
we refer to the microscopic model for Ce-based HF compounds. 

The similarities in the behavior of Ce-based heavy-fermion systems to that of dilute mag- 
netic alloys have led to the assumption that these systems are "Kondo lattices" where the ob- 
served anomalous behavior can be explained in terms of periodically repeated resonant Kondo 
scattering. This ansatz provides a microscopic model for the formation of a singlet ground- 
state and the existence of heavy quasiparticles. An extensive discussion is given in [4]. Direct 
evidence for the Kondo scenario comes from photoelectron spectroscopy. The characteristic 
features of a Kondo system can be summared as follows [14, 15]: At high temperatures, the 
combined PES/BIS spectra from photoemission and inverse photoemission exhibit two distinct 
peaks below and above the Fermi energy. These two features correspond to the valence transi- 
tions f — > f"^\ respectively. The changes in the occupation of the Ce 4f- shells are associated 
with energies of order eV. The high-temperature state can be modelled by weakly correlated 
conduction electrons which are weakly coupled to local f-moments. The f-derived low-energy 
excitations are those of a system of local moments. The direct manifestation of the low-energy 
scale is the appearance of a sharp peak near the Fermi energy at low temperatures. In Ce sys- 
tems, this many-body feature known as "Abrikosov-Suhl" or "Kondo" resonance, is centered 
at Cf+kT* slightly above the Fermi edge Cf. The "Kondo temperature" T* determines the en- 
ergy scale of the dynamical screening of the impurity spin by conduction electrons [4]. The 
evolution of the Kondo resonance with temperature was recently observed by high-resolution 
photoemission experiments [16]. 

The resonance is a genuine many-body feature reflecting the small admixture of f^ con- 
figurations to the ground state and the low-lying excitations which are mainly built from f^- 
configurations. At sufficiently low temperatures T -C T*, the contribution of the narrow reso- 
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nance peak to the thermodynamic and transport properties can be described in terms of a Landau 
theory with heavy fermionic quasiparticles [17]. Based on the corresponding effective Hamil- 
tonian Nozieres [18] introduced a resonant phase shift to account for the impurity contribution 
to the low-energy properties. 

The novel feature observed in stoichiometric Ce-compounds is the formation of narrow 
coherent bands of low-energy excitations. The heavy fermions arise from a decoherence- 
coherence crossover occurring at low temperatures. 

The calculation of realistic quasiparticle bands in Ce-based Heavy Fermion compounds 
proceeds in several steps. For a detailed description see [19]. The first step is a standard LDA 
band structure calculation by means of which the effective single-particle potentials are self- 
consistently generated. The calculation starts, like any other ab-initio calculation, from atomic 
potentials and structure information. In this step, no adjustable parameters are introduced. The 
effective potentials and hence the phase shifts of the conduction states are determined from 
first principles to the same level as in the case of "ordinary" metals.The f-phase shifts at the 
lanthanide and actinide sites, on the other hand, are described by an empirical resonance type 
expression 

T]y~tan"^ — ^— (2) 

charactered by the position ey. and the width of the many -body resonance. One of these 
free parameters is eliminated by imposing the condition that the charge distribution is not sig- 
nificantly altered as compared to the LDA calculation by introducing the renormalation. The 
renormalised band method devised to calculate the quasiparticles in heavy-fermion compounds 
thus is essentially a one-parameter theory. We mention that spin-orbit and CEF splittings can 
be accounted for in a straight- forward manner [20]. 

The strong local correlations in Kondo lattices lead to an observable many-body effect, i. e., 
the change with temperature of the volume of the Fermi surface. At high temperatures, the /- 
degrees of freedom appear as localised magnetic moments, and the Fermi surface contains only 
the itinerant conduction electrons. At low temperatures, however, the / degrees of freedom are 
now tied into itinerant fermionic quasiparticle excitations and accordingly, have to be included 
in the Fermi volume following Luttinger's theorem. Consequently the Fermi surface is strongly 
modified. This scenario [21] was confirmed experimentally by measurements of the de Haas- 
van Alphen (dHvA) effect [22-24]. 

2.2 DUAL MODEL FOR U-BASED SYSTEMS 

The Kondo picture, however, does not apply in the case of the actinide compounds. The dif- 
ficulties with this model have been discussed in [25]. The difference between the Ce-based 
heavy-fermion compounds and their U-counterparts can be seen directly from the photoemis- 
sion spectra [14]. In U-based heavy-fermion compounds, the fingerprint character of the transi- 
tions f is lost. Instead of the well defined f-derived peaks familiar from the Ce systems 
we encounter a rather broad f-derived feature. This fact shows that the f- valence in the actinide 
heavy-fermion systems is not close to integer value as it is the case in Ce-based HF compounds. 
In fact, the f-valence of the U ions has been discussed rather controversially. 
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Increasing experimental and theoretical evidence points towards a dual nature of the 5/- 
electrons in actinide -based intermetallic compounds, with some of the 5/ electrons being lo- 
calised and others delocalised. The dual model which assumes the coexistence of both delo- 
calised and localised 5/-electrons provides a scheme for microscopic calculations of the heavy 
quasiparticle bands in actinide compounds. The concept is summarised as follows: The "de- 
localised" 5/-states hybridise with the conduction states and form energy bands while their 
"localised" atomic-like counterparts form multiplets to reduce the local Coulomb repulsion. 
The two subsystems interact which leads to the mass enhancement of the delocalised quasipar- 
ticles. The situation resembles that in Pr metal where a mass enhancement of the conduction 
electrons by a factor of 5 results from virtual crystal field (CEF) excitations of localised Af^ 
electrons [26]. The underlying hypothesis is supported by a number of experiments. For an 
extensive discussion we refer to [27]. Calculations based on the dual model as described above 
reproduce the dHvA data in UPts [28] and UPd2Al3 [29]. The results of the latter calculation 
will be discussed in sect. 5.2. 

The dual nature of the 5/ states which is found in many actinide intermetallic compounds 
is a consequence of the interplay between local Coulomb correlations and hybridisation with 
the conduction electrons. The underlying microscopic mechanism is an area of active current 
research (see [30] and references therein). LDA calculations show that the effective hopping 
matrix elements for different 5/ orbitals are unequal. But it is of interest to understand why 
only the largest among them is important and why the others are suppressed. 

The reason lies in the competition between anisotropic hybridisation and angular correla- 
tions in the 5f shell which can be seen by exact diagonalisation of small clusters modelling the 
U sites in heavy fermion compounds [30]. Keeping only the degrees of freedom of the 5f shells 
while accounting for the conduction states by (effective) anisotropic 5f-intersite hopping leads 
to the model Hamiltonian 

^ = -1 tj^{cl{n)cj^{m)+h.c)+Hc (3) 

The first sum is over neighboring sites {nm). Furthermore ft (n) (cy,(n)), creates (annihilates) 
an electron at site n in the 5/ j = 5/2 state with = —5/2, ... ,5/2. The effective hopping 
between sites results from the hybridisation of the 5/ states with the orbitals of the ligands and 
depends generally on the crystal structure. Rather than trying to exhaust all possible different 
lattice symmetries, we shall concentrate here on the special case that hopping conserves j^. 
While this is certainly an idealisation, it allows us to concentrate on our main topic, i. e., a study 
of the influence of atomic correlations on the renormalisation of hybridisation matrix elements. 
The parameters tj^{— t^jj are chosen in accordance with density-functional calculations for 
bulk material which use jj^ basis states. The local Coulomb interactions can be written in the 
form 

where the Coulomb matrix elements f/jv,j,27 3>4 expressed in terms of the expectation val- 
ues Uj of the repulsion between electron pairs in states with total angular momentum J. The 
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Figure 1: Left panel: Enhancement of hopping anisotropies due to intra-atomic correlations in 
a two-site cluster. The orbital-projected expectation value of the kinetic energy Tj^/tj^ along the 
lines connecting linearly the points is written below the figure. Regions with =15/2, 5/2 and 
3/2, are labeled I, II and III respectively. Right panel: Dispersion of the low-energy peak in the 
spectral functions (Ay,(k, co) + A_j,(k, to)) for strongly anisotropic systems \t3/2\ = 

\h/2\ — \h/2\) as well as in the isotropic limit (|?3/2| = \t\ — \t'\ — \h/2\ — \h/2\) 
parameter set of hopping parameters, the spectral functions for k=0 and k=7t are compared. The 
dotted curve denotes the isotropic peak in the atomic limit. 
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actual calculations U j values used in the calculation were determined from LDA wavefunc- 
tions for UPts [28], i. e., Uj=4 = 17.21 eV, Uj=2 = 18.28 eV, and Uj=o = 21.00 eV. We expect 
U 7=4 < U 7=2 < U 7=0 always to hold for Coulomb interactions, independently of the chemical 
environment. In contrast, the relative order of the hopping matrix elements will vary strongly 
from one compound to the next. The average Coulomb repulsion is irrelevant for the low-energy 
physics of the model. It simply restricts the relevant configurations to states such that each site 
is occupied either by 2 or 3 f electrons. The low-energy sector is exclusively determined by the 
differences of the Uj values, which are of the order of leV and thus slightly larger than typical 
bare / - band widths. The latter are obtained, e.g., from LDA calculations for metallic uranium 
compounds like UPts. To mimic the situation in the U-based heavy-fermion compounds we 
consider the intermediate valence regime. The Hamiltonian eq. (3) conserves — L„/^(n) 
where is the z-component of the total angular momentum of the system and the Jz{n) re- 
fer to angular momentum projections on individual sites. We shall therefore characterise the 
eigenstates by their value. Strong on-site correlations result in a considerable enhancement 
of anisotropies in the bare hopping matrix elements. This can lead to a localisation of electrons 
in orbitals with relatively weak hybridisation. The latter is effectively reduced to zero in those 



cases. 

The degree of localisation or, alternatively, of the reduction of hopping of a given orbital 
by local correlations, is quantified by the ratio of the j^- projected kinetic energy Tj^ and the 
bare matrix element tj^ 

(^g.|(4.(n)c±,-(m)+/j.c.)|^g.). (5) 

The ground-state wavefunction \^gs) contains the strong on-site correlations. A small ratio of 
Tj^/tj^ indicates partial suppression of hopping for electrons in the ± orbitals. Two kinds of 
correlations may contribute to that process. The first one is based on the reduction of charge 
fluctuations to atomic and configurations. This is a result for large values oi Uj and can 
be studied by setting all U j equal to a value much larger than the different tj^. The second one is 
due to differences in the U j values, i. e. , ^7 j=4 < U j=2 < U j=o. The differences in the U j values 
are the basis of Hund's rules. Hopping counteracts Hund's rule correlations and vice versa. 
What we want to stress is the fact that those correlations can lead to a complete suppression of 
hopping channels except for the dominant one which shows only little influence. 

Results for the ratios Tj^/tj^ are shown in fig. 1 for a two-site model [30]. As the relevant 
correlations are local the general results qualitatively agree with those found for a three-site 
cluster and four-site clusters [31]. We can distinguish three different regimes with = 15/2, 
5/2 and 3/2, labeled I, II and III respectively. One observes that in region I only the dominant 
hybridisation of the = 3/2 orbital survives while that of the = 1/2 and = 5/2 orbitals is 
completely suppressed. On the other hand in regions II and III the correlation effects on different 
orbitals are are not very different. These findings demonstrate that in particular Hund's rule 
correlations strongly enhance anisotropics in the hopping. For a certain range of parameters this 
may result in a complete suppression of the effective hopping except for the largest one, which 
remains almost unaffected. This provides a microscopic justification of partial localisation of 
5/ electrons which is observed in a number of experiments on U compounds and which is the 
basis for further model calculations described later on. 

Many aspects of partial localisation in the ground state are described appropriately by both 
a Gutzwiller type variational wave function and by a treatment which keeps only those interac- 
tions which are present in LDA+U calculations [32]. This finding should encourage further ap- 
plications of LDA+U and related approaches as SIC-LDA for ground-state properties of actinide 
heavy-fermion materials. The subtle angular correlations which determine the magnetic char- 
acter are accounted for only in limiting cases. In addition, we cannot expect these schemes to 
reproduce the small energy scales responsible for the heavy fermion character of the low-energy 
excitations. The dual nature of the 5/ states is manifest in the spectral functions Ay^(k, (o) and 
the density of states (DOS) A^(tt)) given by 

A^(«) = II^A(k,«) . (6) 
k k 

The spectral functions which are defined by 

n 
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Figure 2: Density of states for a linear 
chain calculated within CPT. The short- 
range correlations are treated within a two- 
site cluster. The hopping parameters are 



= 0.2 eV and \t 
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= 0). 



The high-energy side-bands have predomi- 
nantly = ±5/2 character. 



measure the weights for removing an electrons in state (k j^) from the ground state (negative 
energies) or adding on to it (positive energies). Here | Yq ) denotes the exact ground state of 

an AT-particle system with ground state energy £'q^^ while | Yn^^ ) and | V^^^ ) are complete 

sets of eigenstates of the 1- and — 1 -particle systems with energies E^'^^^ and E^~^\ 
respectively. Technically, the spectral functions are obtained from the one-particle Green's 
functions according to 

Aj^ (k, (O) = - -ImGj^^ (k, (o) (8) 

which are approximately calculated for a cluster following [33]. As a first step towards a micro- 
scopic treatment of partial localisation in extended systems the one-particle Green's functions 
are calculated within cluster perturbation theory. The method combines exact diagonalisation 
of small clusters with strong-coupling perturbation theory [34]. The many-particle states of the 
clusters are determined by means of the Jacobi-Davidson method [35]. The formation of the 
heavy quasiparticles, i. e., the enhancement of the effective mass over the bare band mass is 
related to the transfer of spectral weight from low energies to high-energy satellites. Differ- 
ent kinds of processes may contribute to this phenomenon. First, there is the suppression of 
charge fluctuations which effectively restricts the ground state in the mixed- valent regime to 
and configurations. Neglecting angular correlations, the spectral function will have a low- 
energy peak from the transitions and whose weight will be reduced due to 
the atomic valence transitions ^ and f^- These transition appear in the spectrum 
at energies well separated from the Fermi level. Second, there are the angular correlations as 
described by Hund's rules. They transfer spectral weight to atomic side bands which appear at 
energies ~ 1 — 2 eV below (or above) the Fermi edge. These structures can be seen from the 
DOS displayed in fig. 2. Experimentally, the "Hund's rule structures" are observed as "humped 
features" in photoemission experiments [36]. Due to the atomic transitions the overall width is 
much larger than the value predicted by standard band structure calculations. 

Suppression of charge fluctuations and angular correlations reduce the weight of the low- 
energy peak to 0.41. This reduction corresponds to an isotropic enhancement of the quasipar- 
ticle mass by ~ 2.5. The quasiparticle weight is further reduced by hopping between the sites. 
The latter introduces dispersion effects and orbital-dependence as can be seen from fig. 1 (right 
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panel). The spectral weight of the isotropic "atomic" peak is distributed over a rather broad en- 
ergy range. The strong local correlations result in a considerable enhancement of anisotropics 
in the bare hopping matrix elements. The channels with the dominant hopping exhibit well- 
defined dispersive quasiparticle peaks in the low-energy regime whereas the remaining channels 
contribute a non-dispersive incoherent background. The results are consistent with previous re- 
sults the orbital-projected expectation values of the kinetic energy. 



2.3 THE NFL STATE AND ITS SCALING LAWS 

In the vicinity of a second order phase transition the behavior of various physical properties ex- 
hibit singularities which can be characterised by critical exponents [37]. The latter are universal 
in the sense that they do not depend upon the detailed microscopic nature of the system. They 
are determined by the dimensionality of the system and the degrees of freedom associated with 
the long-range correlations in the ordered phase. As a consequence, the critical exponents have 
been successfully used to classify the critical behavior in classical systems. The universality 
reflects the fact that the characteristic length and time scales, i.e., the correlation length ^ and 
the relaxation time T diverge like 

^ - l^r'' and T ~ l^r''" (9) 

for |g| (see e. g. [38, 39]) where g = t, r measures the distance to the critical point according 
to 

T — Tr p — Pc H — He , 

r = — ^, r=^^-^ or r=——^ (10) 

Pc 

Here T, p and H are temperature, pressure and field while Tc, pc and He are their critical 
values. It is important to note that the critical exponents are not independent. They are related 
by scaling relations and, in addition, they are connected with the dimension d of the system 
by hyperscaling relations. On approaching T^ from above (t^ 0+) the growing amplitude of 
thermally excited collective modes drives the system to a new state with spontaneously broken 
symmetry characterised by an order parameter. 

On the other hand quantum phase transitions QPT) which take place at zero temperature 
(Tt =0) are driven by the contribution of quantum fluctuations to the ground state energy which 
depends on a control parameter r ^ 0+ associated with pressure or field. In this case the 
effective dimensionality dg/f is given by 



deff = d + z (11) 

where z is the dynamical exponent introduced above. Due to the effective increase in dimension- 
ality the latter may reach the upper critical dimension. As a consequence, the critical exponents 
associated with the QPT may assume mean-field values [40]. In general, the characterisation of 
a continuous quantum phase transition is more subtle than that of a classical one. 

Besides specific heat, susceptibility and resistivity the Griineisen ratio F has turned out to 
be a sensitive quantity to charaterise the vicinity of the QCP. Zhu et al. [41] have shown that the 
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r ratio defined in terms of the molar specific heat Cp = ^ {dS/ dT)p and the thermal expansion 
a = ^ {dV/dT)p^ = 1 {dS/dp)T^j^, namely 

g _ 1 ids/ dp) J 

Cp VmT{dS/dT)p ^ ^ 

diverges at any QCP. Here S is the entropy and V»j = ^ is the molar volume. Away from the QCP 
when the low-temperature behavior is characterised by a single energy scale E* the Griineisen 
ratio in eq. (12) reduces to 

1 dE* 

r~-— ^ (13) 
VmE* dp 

This form was already introduced and discussed in [42-45]. Instead of pressure, an external 
magnetic field H may be used as control parameter. In this case, the ratio 

^ CH\dT)^ T{dS/dT)^ T\dH)s ^ ^ 

can be determined directly from the magneto-caloric effect. The Griineisen ratio Y and its 
magnetic counterpart Th owe their importance to the fact that at the critical value r = the 
temperature dependence is given by 

r(r,r = o)-r-^ (15) 

and consequently permits to measure vz. Comparison with experiment usually requires to ac- 
count for corrections to the simple scaling ansatz from which eq. (15) is derived. The correc- 
tions, however, have to be evaluated from a microscopic model. For a discussion of the validity 
of eq.(15) we refer to [46, 47]. 

Neutron diffraction results for the A-phase in CeCu2Si2 suggest that the magnetic quantum 
phase transition in stoichiometric heavy fermion compounds arises from the instability of the 
strongly renormalised Fermi liquid with respect to the formation of a SDW. Model studies for 
the SDW transition therefore deserve special attention. Using the Ginzburg-Landau-Wilson 
functional for the effective action [40] 

5[0] = I ^+?'+^')|^q.-a^|'+5(^^ with rq = ro^-2 

5(4) = u J d'^r\^(r,T)f (16) 

the thermal expansion and Griineisen ratio were calculated on the nonmagnetic side of the 
phase diagram [41]. Here 4> is the order parameter that fluctuates in space and imaginary time, 
/3 = I /kT is the inverse temperature and (On — 2nnT are bosonic Matsubara frequencies. The 
parameter r controls the distance from the critical point. The Landau damping which is linear 
in Icq I is due to the scattering of quasiparticles by spin fluctuations. It is characterised by the 
dynamical exponent z. The theory starts on the Fermi liquid side (r > r^) which is characterised 
by a large Fermi surface including the f-degrees of freedom. The low-energy excitations are 
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Figure 3: Thermal expansion showing the 
supression of NFL behaviour as function of 
field. The inset shows that for B = the 
critical contribution to the Griineisen ra- 
tio of CeNi2Ge2 (sect. 3.2) scales like T ~ 
1/T-* with X = 1. After eq.(15) this means 
(assuming z = 2 for AF SDW) a mean field 
correlation length exponent v = 1/2 which 
is in agreement with dg// = d+z = 5 for the 
effective dimension [48]. 

1 2 3 4 5 6 
T(K) 

fermionic quasiparticles and their collective excitations. Close to the QCP the static suscepti- 
bility is assumed to diverge at a specific wave vector Q which signifies the transition into the 
SDW state. The quasiparticles are strongly scattered along "hot lines" which are connected by 
Q [49]. This strong scattering modifies the low-temperature thermodynamic and transport prop- 
erties which exhibit anomalous scaling relations close to QCP. They differ from those familiar 
from Fermi liquids. The scaling relations are derived adopting the renormalisation scheme of 
Millis [50, 51] which proceeds in close analogy to the spin-fluctuation theory of ferromagnetism 
[52]. The scaling behavior obtained for various models and many non-Fermi liquid compounds 
are summarised in [53]. 

The model calculations in [41] classify transitions according to various types of magnetic 
order. Two-dimensional (d=2) as well as three-dimensional (d=3) systems are considered as- 
suming the values z = 2, 3 for the dynamical exponents. Comparison of calculated Griineisen 
ratios with experiments in CeNi2Ge2 [48] shows that the QCP in this stoichiometric HF com- 
pound are consistent with the SDW scenario for three-dimensional critical AF spin fluctuations 
(fig. 3). 

3 CERIUM-BASED HF SUPERCONDUCTORS 

Since the discovery of heavy electrons and their superconducting state [54] in CeCu2Si2 the 
Ce-based tetragonal 1 12-compounds (space group I4/mmm, see fig. 4) have been at the focus of 
investigation. This class of materials offers a unique stage to watch the interplay of quantum- 
criticality, heavy Landau Fermi liquid (LFL) and non- Fermi liquid (NFL) behaviour as well 
as unconventional superconductivity. Although it has a long history, major aspects of the par- 
ent compound CeCu2Si2 of this class where only understood very recently: Firstly the long 
mysterious ' A-phase' associated with the low pressure quantum critical point, has finally been 
identified as an incommensurate (IC) SDW phase caused by the nesting of heavy quasiparticle 
FS sheets [55]. In fact it is now possible to follow the continuous evolution of the SDW wave 
vector and moment size of magnetic phases in CeCu2(Sii_;cGe;,^)2 from the small moment IC- 
SDW in CeCu2Si2 to the atomic moment commensurate antiferromagnetism in CeCu2Ge2 [56]. 
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Secondly the anomalous stability of SC in CeCu2Si2 for pressures far above the A-phase 
QCP [57] which did not exhibit the typical dome-shape around the QCP pressure has now been 
understood. Likewise CeCu2Ge2 does not have a dome-shaped SC Tc(p) curve once AF order is 
suppressed. Employing a negative pressure shift caused by the larger Ge- radius it looks rather 
similar to that of CeCu2Si2. Artificial reduction of the SC condensation energy in CeCu2Si2 by 
Ge- alloying and applying positive shifts with hydrostatic pressure leads to a breakup of the SC 
region into two isolated domes which are associated with the previous A-phase magnetic QCP 
and a new high pressure valence transition QCP respectively [58]. This also suggests that SC 
phases and pairing mechanisms in the isolated SC regimes are different. Quantum criticality in 
the CeCu2(Sii_xGex)2 system is connected with distinct NFL anomalies that can only partly 
be explained within the standard nearly AF Fermi liquid picture, notably the discrepancy of the 
temperature scaling of resistivity and 7- coefficients at higher fields close to the A- phase QCP 
may require different concepts. 

The HF compounds CeRh2Si2 and CePd2Si2 also become superconducting only under ap- 
plication of pressure like CeCu2Ge2. As in CeCu2Ge2 AF order in CeRh2Si2 vanishes in a first 
order phase transition at the critical pressure p^. SC in the latter appears only in an extremely 
narrow region around p^, unlike in the former compound. Due to the first order transition both 
CeCu2Ge2 and CeRh2Si2 are in a LFL state close to p^ For CePd2Si2 on the other hand a 
QCP where AF order vanishes continuously leads to a SC p-T phase diagram which is of the 
canonical single dome shape centered at pc where clear cut NFL behaviour in the resistivity 
is observed. Finally CeNi2Ge2 is quite unique because it is not magnetic but rather displays 
almost ideal NFL behaviour at ambient pressure although incipient superconductivity is found 
at very low temperatures in highest-purity samples. 

A comparison of the various Cel22 compounds leads one to conclude that the existence of 
an AF QCP with associated NFL behaviour is sufficient to induce superconductivity around pc. 
Indeed the AF spin fluctuation theory (sect. 4) predicts a stabilisation of the SC state when the 
QCP is approached. However, as the examples of CeCu2Ge2 and CeRh2Si2 show, SC in Cel22 
compounds may just as well appear in the normal LFL state. 

It took more than 20 years to discover an additional class of Ce-based HF superconductors 
with the general formula Ce„MmIn3„_|_2m (M = Co, Ir and Rh). This belongs to the tetragonal 
layer type materials derived from the 'infinite layer' parent compound Celn3 by stacking n 
'Celns' and m 'MIn2' subunits along the tetragonal axis. While CeCoIns and Celrlns are single 
layer (n=l) HF compounds which are superconducting at ambient pressure [59], CeRhlns is 
a local moment AF which becomes SC only under pressure. Likewise in the bilayer (n=2) 
compounds SC appears both at ambient pressure for Ce2CoIn8 and Ce2lrln8 and at finite pressue 
for Ce2RhIn8. Especially CeCoIns has turned out to be of great interest. Its comparatively large 
Tc has allowed to determine the d-wave symmetry of the SC gap function by field-angle resolved 
magnetothermal conductivity [60] and specific heat [61] measurements. High-field specific heat 
[62] and ultrasonic attenuation measurements [63] have indicated that firstly the SC transition 
evolves to first order at larger fields, and secondly evidence for the long-sought FFLO phase 
with SC pairing at finite momentum exists. As in the Cel22 class SC appears preferably in the 
vicinity of magnetic QCPs and the corresponding AF/SC phase diagrams have been determined 
in pressure experiments. Likewise NFL phenomena above the AF QCP have been observed. 

Another exotic HF Ce- superconductor, CePtsSi has now been discovered which is the first 
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Figure 4: Left: Tetragonal unit cells 

of the ThCr2Si2 and CeM2X2 (Cel22) 
structure where M = Cu, Ni, Ru, Rh, 
Pd, Au, .. (black circles); X = Si, 
Ge (small grey circles) and Ce (large 
grey circles). Right: Conventional 
unit cell of the CeMlns (Cel 15) struc- 
ture whith M = Co, Ir, Rh (black cir- 
cles); In (small grey circles) and Ce 
(large grey circles). 
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with a non-centrosymmetric structure [64]. Superconductivity in crystals without inversion cen- 
ter was discussed theoretically first in [65]. The lack of inversion symmetry raises fundamental 
questions of classification of the SC states according to odd and even parity [66]. It remains to 
be seen whether this compound is the first member of a new class of HF SC. 



3. 1 QUANTUM CRITICAL BEHAVIOUR AND SUPERCONDUCTIVITY IN CeCu2Si2AND 
THE ALLOY SERIES CeCu2(Sii_;,Ge;,)2 

The first Ce- based HF metal and superconductor [54] has for a considerable time served as 
a most fruitful model for investigating strongly correlated electron systems. Until rather re- 
cently it was also the only HF Ce compound exhibiting superconductivity. Most importantly 
in CeCu2Si2 Si can be continuously subsituted by Ge while the compound changes from a HF 
superconductor to a normal antiferromagnet [56]. The Ge substitution acts like a negative chem- 
ical pressure, likewise fully substituted CeCu2Ge2 under positive hydrostatic pressure behaves 
similar as CeCu2Si2 under ambient pressure. Ge- substitution and application of hydrostatic 
pressure therefore allows to vary the coexistence and competition behaviour of superconducting 
and other order parameters in the alloy series CeCu2(Sii_;<:Gex)2 • 

Indeed already stoichiometric CeCu2Si2 at ambient pressure shows an additonal 'A-phase' 
which envelops SC in the B-T phase diagram [67] (fig. 5). When T^- and T^ are close they do 
not coexist and the A-phase is expelled from the SC region of the B-T plane. Ge-substition 
or suitable deficit of Cu however stabilises the A-phase (T^ > T^) and coexistence of both 
phases results. Hydrostatic pressure and Cu excess on the other hand reduces T^ to zero at 
the quantum critical point (QCP). Even after single crystals became available [68] the nature 
of the A-phase has remained mysterious. Later, neutron diffraction on weakly Ge- substituted 
CeCu2(Sii_;(:Gex)2 samples have shown that the magnetic order observed in the stoichiomet- 
ric CeCu2Ge2 compound is still present, although with continuously decreasing moment [69]. 
The slightly incommensurate (IC) modulation vector did not change very much as function of 
X (Ge concentration). Finally IC magnetic order with Q = (0.215, 0.215, 0.530) (r.l.u) has 
been found in the stoichiometric compound CeCu2Si2 by neutron diffraction [55], which iden- 
tified the A-phase as a conventional IC-SDW phase with rather low moment ~ 0. 1 jib per Ce 
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7 [mJ/molK^] 


Ta,[K] 




Tc.(p=0)[K] 


Tr(Pm)[K] 


CeCu2Si2 


700 


0.8 


0.10 


0.7 


2.3 (3 GPa) 


CeCu2Ge2 


1200(1) 


4.3 


1.05 




1.7 (16 GPa) 


CeNi2Ge2 


350 








0.4 (2.2 GPa) 



Table 1: Material parameters of Ce- based HF superconductors. Value extrapolated to T = 
by fitting the S=l/2 Bethe ansatz result to data for T > T;v using T* = 6K. 

site (fig. 6). The ordering wave vector nicely agrees with the nesting wave vector of heavy FS 
sheets in CeCu2Si2 found by Zwicknagl et al. [20, 55]. It is clear now that the A-phase of 
CeCu2Si2 evolves continuously from the magnetic phase of CeCu2Ge2 , when x is reduced 
from X = 1 progressively . Taken together the pressure and Ge- substitution dependences at 
first sight do not support the idea that SC appears around the QCP of the A-phase because the 
Tc(x,p) curve (full dots in fig. 7) lacks the typical 'dome shape'. Specifically Tdp) strongly 
increases for increasing hydrostatic pressure p [57] and reaches its maximum far away from the 
QCP (fig. 7). This feature is very different from the conventional SC/AF QCP phase diagram 
which applies for other Ce- based HF superconductors like CePd2Si2 and CeRh2Si2 (fig. 16). 
However it is now clear [58] that this anomalous high pressure maximum in T^ is due to the 
presence of an additional high pressure QCP associated with an only weakly first order valence 
phase transition of Ce which does not break any symmetry [70]. The charge fluctuations asso- 
ciated with the valence transition have been proposed to mediate the SC pairing around the high 
pressure maximum [71] which is distinct from both spin fluctuation and magnetic exciton mech- 
anism invoked for other Ce- and U-HF compounds respectively. The appropriate dome shape 
of two separated QCPs has been revealed by using simultaneous tuning with Ge-substitution 
and hydrostatic pressure. The latter compensates the effect of the former, but the mean free 
path is decreased which destabilises the unconventional SC order parameter and reduces it's 
T^ . In this way the large stability region of stoichiometric CeCu2Si2 may be broken into the 
two domes associated with the magnetic and valence transition QCPs (fig.7). It seems now 
that CeCu2Si2 and substitutes also fit into the generic QCP scenario of Ce-HF superconductors 
albeit with the surprise of an additional qualitatively new QCP. 

3.1.1 QUANTUM CRITICAL POINTS AND THE B-T PHASE DIAGRAM 

Like the other Cel22 compounds the CeCu2(Sii_;cGe;c)2 series crystallises in the tetragonal 
body centered ThCr2Si2 structure (fig. 4). The stability of A and SC phases is extremely sensi- 
tive to the stoichiometry of the constituent elements. The homogeneity range of CeCu2Si2 in the 
ternary chemical Ce-Cu-Si phase diagram is quite limited and comprises three regions where ei- 
ther A or SC alone are present or A/SC coexist [72]. This agrees with the observation that in the 
combined Ge- doping and pressure phase diagram shown in (fig. 5) the crossing of T^ and T^ 
curves is located around the stoichiometric case at ambient pressure. Around the first low pres- 
sure magnetic QCP where Ta ^ NFL behaviour with Ap{T) ~ T^/^ and y{T) = aT^/^ 
was observed for moderately low temperatures and moderately high fields [75]. This low tem- 
perature phase diagram which has been composed of different transport and susceptibility mea- 
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Figure 5: Left panel: Coexistence/competition phase diagram of SC and magnetic A phases 
close to the magnetic QCPl. Phase boundaries are obtainded from susceptibility and resistivity 
measurements and composed of alloying (negative pressure) and hydrostatic pressure results, 
separated by a dashed line. Here g ~ 1-x on the left and g ~ p on the right [72] (c.f. theoretical 
results in fig. 10). Right panel: B-T phase diagram of CeCu2Si2 for B || a. Original version in 
[73], completed version from [74]. For this sample the A-phase is expelled from the SC region 
(no coexistence). 



surements already shows that Tc saturates in a plateau after the T^- crossing instead of forming 
the canonical QCP dome shape. Indeed on increasing pressure even further Tf increases steeply 
reaching a maximum around 3 GPa [76] to 4 GPa [57] (fig. 7). The delicate competition and 
coexistence behaviour for stoichiometric single crystals of CeCu2Si2 has already been obvi- 
ous from earlier investigations on the B-T phase diagram obtained with ultrasonic methods and 
dilatometry [73, 74] which is shown in fig. 5 (right panel). For this sample Tc < T^ and A/SC 
phases do not coexist. The expulsion of the A-phase from the SC region in fig. 5 can be clearly 
seen from sound velocity anomalies across the A and SC phase boundaries. The still poorly 
understood B phase might be another SDW phase with different propagation vector Q and/or 
SDW moment polarisation. 



Figure 6: Neutron diffraction intensity in 
reciprocal (hhl) plane in CeCu2Si2 at tem- 
peratures above and below the A-phase 
transition temperature T^. The magnetic 
peak appears at (0.215,0.215,1.470) cor- 
responding to an incommensurate SDW 
modulation vector Q = (0.215,0.215,0.530) 
[55]. 




16 



Figure 7: p-T phase diagram of 
CeCu2(Sii_xGeY)2 composed from 
pressure (p) dependence and Ge-doping 
(x) dependence of p(T), C(T) and xO^)- 
For different substitution level x, the 
pressure is shifted by the value of its 
corresponding critical pressure PcKx) for 
the magnetic QCPl (Ap = p - pci(x)). Two 
SC domes are associated with QCPl(pci) 
(left dome) of the A- phase and QCP2(pc2) 
of the integer valence (Kondo) phase (right 
dome) [58]. The dashed line indicates the 
valence transition in CeCu2Ge2 and T^ 
(dots [76]) corresponds to stoichiometric 
CeCu2Si2. 

Ap (GPa) 

To unravel the puzzling anomalous stability of T^ at high pressure systematic studies with 
slightly Ge doped single crystal CeCu2Si2 were performed [58]. It has been long suspected that 
the anomalous Tc(p) behaviour is associated with the presence of a second QCP connected with 
a Ce-valence transition [70]. To support this scenario one has to prove that the region around 
Tc — Ta and around the maximum Tc correspond to distinct SC phases. This is achieved by 
the following strategy [58]: i) Due to larger Ge ionic radius slight Ge doping corresponds 
to application of internal negative pressure which may be compensated again by application 
of external positive pressure, returning to zero effective pressure, ii) Because the electronic 
mean free path has decreased due to scattering from Ge impurities, the condensation energy 
and hence Tc of the anisotropic SC phase will be decreased, destabilising SC. In this way 
the continuous SC region of stoichiometric CeCu2Si2 indeed breaks up into two separate SC 
domes as shown in fig. 7. The first one centered at Pci(x) is the magnetic QCP associated 
with the A-phase and its critical line Tyv(Ap), the second QCP at Pc2(x) is supposed to describe 
the Ce^+ZCe^^ valence transition where the 4f electron of Ce^+ becomes delocalised. This 
transition does not break any symmetry and therefore its critical line Tv(Ap) (not shown in 
fig. 7) is assumed to end in a critical endpoint. The critical line and its endpoint, however, have 
sofar not been confirmed for these doped CeCu2(Sii_;cGex)2 systems by independent methods. 
However for pure CeCu2Ge2 a weak, symmetry-conserving valence transition was observed via 
x-ray diffractometry [77] nearp ~ 15 GPa, where also Tc(p) assumes it's maximum value [57]. 
Certainly a large volume collapse signifying the 4f delocalisation as in the a — 7 Ce transition 
would not be compatible with the presence of the SC phase around Pc2(x). 

3. 1 .2 THE A-PHASE AND ELECTRONIC STRUCTURE OF CeCu2Si2 

The incommensurate magnetic order of the A-phase suggests a SDW origin associated with 
the itinerant heavy electron FS in CeCu2Si2. To support this conjecture the static magnetic 
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susceptibility xi^l) has to be calculated within the context of renormalised band theory [20]. 
In this approach (sect. 2) resonant phase shifts r}j>^ for the 4f-states are introduced empirically 
into the LDA scheme to generate the heavy bands at the Fermi level while charge neutrality 
and proper f-count are preserved. For the three CEF split Kramers doublets of Ce^+ (J=§) the 
ansatz for the phase shift is 

fm 

Here e^^^^ defines the centers of heavy bands with denoting the renormalised 4f- level and 
A,n the energies of the three Kramers doublets taken from neutron scattering. The resonance 
width Fy^ is an empirical parameter which determines the width kT* or mass m* of the heavy 
electron band. It is adjusted to reproduce the proper experimental value of 7 = C/T. The main 
heavy electron sheet obtained in [19] is shown in fig. 8. It consists of stacked columns along c. 
Obviously there are large flat parts on the columns with a nesting vector Q as indicated. Indeed 
this leads to a pronounced maximum of ;f (q) at q = Q as shown in the contour plot of fig. 8. The 
experimental IC ordering vector of the A-phase (fig. 6) agrees well with the maximum position 
of X (q) giving credence to the SDW picture for the A phase. 

1 n 




Figure 8: Left panel: Main heavy FS sheet in CeCu2Si2 where columnar nesting with wave 
vector Q is indicated. Parameters for renormalised band calculations are: T* ~ 10 K and A(-igp 
= 330 K overall CEF splitting. Right panel: Comparison of experimental propagation vector 
(fig. 6) and contour map of theoretical xin) in (hhl) plane [55]. 



3.1.3 THE HIGH PRESSURE MIXED VALENT PHASE TRANSITION 

The present experimental evidence for the second QCP at high pressure connected with Ce- 
valence change was described in [70]. A valence transition does not break any spatial symmetry 
but charge fluctuations and volume strain are strongly coupled which commonly leads to a first 
order valence transition as in the famed y-a transition of Ce. No such volume change has been 
found at p^i of CeCu2Si2 making it a rather exceptional case. On the other hand this feature may 
provide a link to superconductivity as discussed in the next section. The only direct evidence for 
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Figure 9: Behaviour of y-coefficient, resid- 
ual resistivity po and A-coefficient in p(T) 
= po+AT^ against ~ T* which scales 
monotonously with pressure [70]. 



a small valence change was obtained earlier in Ljjj- x-ray absorption experiments [78]. There is 
however considerable indirect support for the valence change. The gradual delocalisation of the 
4f^- electron under pressure leads to a less strongly correlated electron state, i.e. the effective 
mass will be reduced according to 

^ = (18) 

m 1 ~ 

when the f-level occupation Wf h reduced significantly below its value in the Kondo limit 
ny^ <1. The linear specific heat coefficient y - C(T)/T should reflect this decrease. Likewise 
the Kadowaki- Woods ratio A/y^ should decrease accross the valence transition [70]. If it is not 
of first order, critical valence fluctuations should also lead to a strong increase in the residual 
resistivity around the transition. These expected fundamental features of a valence transition 
are indeed present in CeCu2Si2 as shown in fig. 9. Theoretical models for the present case of 
a rapid but continuous valence transition under pressure are based on the extended periodic An- 
derson model (PAM) including a Coulomb repulsion U/c between conduction and f- electrons, 
in addition to the on-site U// Coulomb term in the common PAM. Within the more simple 
impurity model this term has the tendency to decrease the f-occupation and when the energy 
of the 4f^ state, ef+\Jfc approaches ep, i.e. the energy of the 4f^ hole state, Wf is rapidly re- 
duced to values much less than one (fig. 1 1). External pressure is thought to tune primarily Uyc 
which is an inter-site Coulomb integral and thus cause the valence transitions. For the lattice 
model the valence change as function of U/c has been calculated within slave boson fluctuation 
approximation [71]. 

3. 1 .4 UNCONVENTIONAL SUPERCONDUCTIVITY AND THE TWO QCP' S 

The breakup of the SC region into two independent domes in fig. 7 allows one to speculate 
that the SC pairing mechanism at each QCP involves exchange of fluctuations of the associ- 
ated order parameters, magnetic (A-phase transition) or charge fluctuations (valence transition) 
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Figure 10: Theoretical phase diagram for t[k] 
A and dj^y- SC order from a two band j 
model. Full and dashed line correspond to 
second and first order transitions. A+SC is 
the coexistence region. The effective tight 
binding hopping integral teff scales with p 
or 1-x (c.f. experimental results in fig. 5) 
[72] 
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respectively. This may possibly imply a different type of superconducting gap symmetry, al- 
though not much is known about the SC state in both regimes, especially for the high pressure 
SC region. In the magnetic QCP regime Cu NQR experiments under pressure for polycrys- 
talline CeCu2(Sii_;cGex)2 [79] and under ambient pressure for Cu and Ce off- stoichiometric 
poly crystalline samples [80] show i) 1/Ti ~ T^ behaviour indicating a SC state with gap nodes 
and ii) the presence of soft magnetic fluctuations close to a pressure or doping regime where the 
SC is stabilised. This suggests that the magnetic QCP scenario for unconventional HF-SC in 
CeCu2Si2 may be correct. Further knowledge about the nature of the SC gap function can only 
come from field-angle resolved pure single crystal investigations, e.g. C(H,T) measurements at 
very low temperatures. 

Model for SC at the magnetic A- phase QCP: 

A 2D two-band model for SDW-SC states in CeCu2(Sii_xGe;c)2 has been studied in some detail 
[72]. The microscopic origin of the pair potential is not specified in the BCS type model. A 
tight binding (TB) band centered at the r(0,0) point and with Q = (5,5) nesting feature repre- 
sents a simplified version of the FS columns in fig. 8 whose nesting properties are known to be 
responsible for the SDW A-phase. In addition spherical electron pockets which represent the 
remaining FS sheets of CeCu2Si2 are included. The former carries both SDW and SC order 
parameters, the latter only SC. Therefore a dxy SC state is favored. The A/SC competition- 
coexistence phase diagram resulting from the coupled gap equations is shown in fig. 10. There 
teff is the hopping element of the TB band which is thought to scale linearly with the control 
parameter g which is 1-x (x = Ge doping) or the pressure p. The comparison of the theoretical 
(fig. 10) and experimental phase diagram (fig. 5) shows a striking resemblance. In this model 
the tail of the coexistence A+SC regime is stabilised by the SC pairing on the spherical FS 
sheets. A more elaborate discussion of this issue using the real FS of CeCu2Si2 is given in 
sect. 3.1.6. 

Model for SC at the valence transition QCP: 

A more microscopic model for the SC mechanism around the high pressure valence transition 
has been proposed in [71]. The extended PAM mentioned above leads to an electron pairing 
mainly via the exchange of c^f charge fluctuations with small wavelength. In weak coupling 
approximation a d-wave state is found to be stabilised in the regime where appreciable charge 
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Figure 11: Valence (n/) transition in the 
extended (impurity) Anderson model as 
function of the f-level position Ef (D = 
conduction band width, = 0). Associ- 
ated Tc dependence of d-wave SC is also 
shown. [71] 



fluctuations, controled by the pressure dependent \Jfc- parameter, set in but the quasiparticle 
DOS (the mass enhancement) is still sizable. The maximum in scales with the steepness 
of the continuous valence transition as seen in fig. 11. This does not incorporate realistic FS 
features but assumes a single spherical conduction band. Also the real change in n/ across the 
valence transition will certainly be much less than in the theoretical model (fig. 1 1). Therefore 
it should be taken as an illustration of the mechanism and not literally as explanation for high 
pressure SC in CeCu2(Sii_;cGex)2. 

3.1.5 NON-FERMI LIQUID STATE AND QUANTUM CRITICAL BEHAVIOUR 

Above the SC domes the two QCP's in CeCu2(Sii_xGe;i;)2 lead to pronounced NFL anomalies 
in thermodynamic and transport properties. Naturally one would invoke here the quasiparticle 
scattering from low energy fluctuations of the associated order parameters [49] as one possible 
mechanism behind NFL behaviour. For the A-phase QCP a detailed comparison of experi- 
mentally observed T,B-scaling laws for p(T) = po + Ap(T), xiT) and /(T) for single crystal 
CeCu2Si2 and polycrystalline samples of various Ge-dopings has been performed [75, 81]. It 
was found that for samples slightly above Ta(x) = the expected NFL scaling laws at a 3D AF 
QCP 

Ap(r)=j8ri Y{T) = ro-ocT-2 (19) 

are indeed well fulfilled. However there are also distinct deviations, notably the high-field 
behaviour of Ap(T) where a crossover of the resistivity exponent from 3/2 to the FL exponent 
2 is observed. In contrast 7(T) remains qualitatively unchanged [75] at moderately low T and 
exhibits an upturn at the lowest temperatures which is not yet understood. Such disparities in 
the scaling of Ap(T) and 7(T) are now also known from other NFL systems and proposals how 
to account for them have been made [75]. The study of NFL behaviour at the high pressure 
valence transition QCP on the other hand is still in its beginning. The NFL behaviour in Ap(T) 
above the two distinct SC domes of fig. 7 is illustrated in fig. 12. 

3. 1 .6 COEXISTENCE OF SUPERCONDUCTIVITY AND A-PHASE IN CeCu2Si2 

The central goal in this section is to study the competition-coexistence phase diagram within 
a more realistic model than above. We start from the real FS sheets of CeCu2Si2 as obtained 
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Figure 12: Left panel: NFL behaviour of Ap(T) in CeCu2(Sii_^Ge;c)2 close to QCPl (left) and 
QCP2 (right) [58]. Right panel: Variation of the SC gap function amplitude |A(k) | / |Ao(r) | = 
|(^r3(k)| for the (pseudo-) singlet wave function with -symmetry <^3(k) ~ cos kxa — cos kya 
on the main sheet of the paramagnetic FS [82, 83]. The amplitude of this SC order parameter 
is maximal on the kidney-shaped surfaces centered along the L-direction which are (almost) 
unaffected by the formation of the A-phase (dark grey). The dominant contributions to the 
latter come from the nesting parts on the heavy columns where the superconducting amplitude 
is small (light gray). 



by renormalised band structure calculations. The FS nesting properties determine the observed 
magnetic structure. Of particular interest are the superconducting states which can either coexist 
with the A-phase or expel it. The model Hamiltonian is given by 

H = Y,ekclci,s+Hint (20) 

ks 

where the first term describes heavy quasiparticles and H,„f their residual interactions. The 
creation (annihilation) operators are for quasiparticles with wavevector k, (pseudo) spins 5 =t 
, I and energy are denoted by cj^^Xck,). The energies which are measured relative to the 
Fermi level are calculated within the renormalised band scheme. The residual interactions in 
the strongly renormalises Fermi liquid is assumed to be repulsive for short separations while 
being attractive for two quasiparticles of opposite momenta on neighboring sites. The former 
favors the formation of a SDW while the latter gives rise to the superconducting instability. 
Within mean-field approximation one has Hint = HsDW+Usc where 

HsDW = -l^«^l^[^(Q7)ick+Q^., + /i.c.] (21) 
Hsc = ^i:[A.,Kk)4Ay+^-^-] (22) 

kss' 

The periodically modulated magnetisation h (Q^) associated with the eight equivalent SDW 
propagation vectors Qj defined below as well as the superconducting pair potential A have to 
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be deteraiined selfconsistently 



A,y(k) = ■^J^giy;yv"(k,k') ^c_,^/^„c,^/^,„^ (24) 

where the strength U of the local Hubbard-type repulsion is of the order of the quasiparticle 
band width ksT* and gss';s"s'" {K k') is the effective pair attraction. The k-summation runs over 
the entire paramagnetic Brillouin zone and L denotes the number of lattice sites. The expec- 
tation values denoted by (...) have to be evaluated with the eigenstates of the total mean-field 
Hamiltonian Hmf — Hq +Hsdw +Hsc, and, consequently, depend upon the order parameters. 
The selfconsistency equations are therefore coupled. The mean-field Hamiltonian implicitly as- 
sumes that the amplitudes of both order parameters are small. In particular, we neglect here the 

pairing amplitudes of the form (^c_^^„c^^q ^,„^. The latter are important when the gaps intro- 
duced by the antiferromagnetic order into the quasiparticle spectrum are large on the scale set 
by superconductivity. For a discussion of this point we refer to [84, 85]. The periodically mod- 
ulated magnetisation associated with the SDW acts on the conduction electrons like a periodic 
spin-dependent potential which we approximate by 

h{x) = Y^h{Qj)e^^r'' (25) 

with identical amplitudes h (Qy) = ho for the eight commensurate wave vectors G {( 

±^ )} which are used for a commensurate approximation to the experimentally found IC 
SDW. The magnetic superstructure breaks the translational invariance of the underlying lattice 
but it conserves the point group symmetry. However the mean-field Hamiltonian is invariant 
under translations with the lattice vectors 

a'l = {2a, 2a, 0) ; 83 = {2a, -2a, 0) ; a'g = {2a, 0, c) (26) 

The volume of the magnetic supercell is 32 times the volume of the paramagnetic unit cell. As 
a result, the Brillouin zone is reduced and the quasiparticle states are modified by extra Bragg 
planes. The opening of new gaps is important at sufficiently low temperatures T <^Ta where 
Ta is the SDW transition temperature. 

Before proceeding with the SC/SDW coexistence analysis of CeCu2Si2 we have to give a 
brief description of the symmetry classification of unconventional SC gap functions. The fun- 
damental property of A^^/(k) is its behaviour as a two-fermion wave function in many respects. 
This expresses the fact that an ODLRO order parameter is not the thermal expectation value of 
a physical observable but rather a complex pseudo- wave function describing quantum phase 
correlations on the macroscopic scale of the SC coherence length. Its phase is a direct signature 
of the broken gauge invariance in the SC condensate. The Pauli principle requires A^y(k) to be 
antisymmetric under the interchange of particles 

A,,,(k) = -A,,,(-k) (27) 



23 



In addition, it transforms like a two-fermion wave function under rotations in position and 
spin space and under gauge transformations. The transformation properties yield a general 
classification scheme for the superconducting order parameter which is represented by a 2 x 2- 
matrix in (pseudo-) spin space. It can be decomposed into an antisymmetric (s) and a symmetric 
(t) contribution according to inA(k) = inAj(k) +mAf (k) with 

3 

mAj(k) = 0(k)im(72 mA,(k) = ^ ii^(k)m(7^/m(72 (28) 

where mCju denote the Pauli matrices. Antisymmetry mA(k) = — inA^(— k) requires 

0(k) = 0(-k) and d^(k} = -d^i-k) (29) 

for the complex orbital functions (^(k) and ^^^(k) (ji = 1-3). For brevity we will frequently 
write A(k) for 0(k) or |d(k)|. 

The order parameter can be chosen either as purely antisymmetric (mAy) or purely symmet- 
ric (mAf) when spin-orbit interaction can be neglected. In the 4f- and 5f- based heavy fermion 
superconductors spin-orbit interaction is strong. As a consequence classification according to 
physical pair spins cannot be used. If their high-temperature crystal structures, however, have 
an inversion center classification according to parity is still possible as is the case in CeCu2Si2. 

The general classification scheme for superconducting order parameters proceeds from the 
behavior under the transformations of the symmetry group ^ of the Hamiltonian. It consists 
of the crystal point group G, the spin rotation group SU(2), the time-reversal symmetry group 
J^, and the gauge group U(l). The appropriate choice of rotations corresponding to weak or 
strong spin- orbit coupling case is determined by microscopic considerations. Using the above 
transformation properties the singlet and triplet gap functions 0(k) and d(k) respectively may 
further be decomposed into basis functions 0p(k) or dp(k) of the irreducible representations 
r (degeneracy index n) of GxSU(2) (weak s.o. coupling) or G (strong s.o. coupling). In the 
latter case the pseudo spin associated with Kramers degeneracy replaces the physical conduction 
electron spin. 

The occurrence of long-range order at a phase transition described by an order parameter is 
most frequently associated with spontaneous symmetry breaking. The simplest superconduc- 
tors where only gauge symmetry is broken are called conventional. In this case the SC order 
parameter has the same spatial symmetry as the underlying crystal, i.e. it transforms as a fully 
symmetric even parity singlet Fi representation of G. It should be noted, however, that conven- 
tional is not a synonym for isotropic, for any G one can form Fi representations from angular 
momentum orbitals of higher order /, for example / > 2 for tetragonal and hexagonal symmetry 
and / > 4 for cubic symmetry. On the other hand, a superconductor with additional broken sym- 
metries besides gauge symmetry is called unconventional. It can have either parity. A recent 
summary is found in the monography [7]. 

NMR experiments [79, 80] in slightly Ge-doped or Cu-off- stoichiometric samples suggest 
that the SC gap is unconventional and has line nodes. In our discussion we restrict ourselves 
to a singlet pair state as revealed by strong Pauli limiting seen in early Hc2 studies [86]. The 
singlet state is characterised by a scalar order parameter 

A(k),,, = 0(k)(/a2),,, (30) 
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The anisotropic effective pair interaction in the (pseudo-) spin singlet channel can be expanded 
according to 



1 

g (k., -k/; -k's", k's'") {iG2),, {i02),,,„ X (rO-)) (k) 0i^> (k') (31) 

j ^ ^ ^ K=\ 

where (k) is a basis function which belongs to the Jcth row, Jc = \,. ..,d^^\ of the d'^^^- 
dimensional representation Y^^^ of the symmetry group. In principle we have to classify the 
order parameters with respect to the antiferromagnetic lattice. Since the two ordering tempera- 
tures are so close we focus on pair states which are compatible with the translational symmetry 
of the paramagnetic lattice. The corresponding basis functions are listed in [5, 19, 87, 88]. In 
the explicit calculations, we restrict ourselves to one-dimensional representations for simplicity. 
The generalisation to multi-dimensional representations is rather straightforward [82]. 

We further simplify the problem by focussing on the symmetry Y which yields the strongest 
quasiparticle attraction g^. Assuming that this most stable order parameter is non-degenerate 
yields the separable interaction kernel 

g (k5, -k/; -k'/', k'/") ^ (/C72),„,,„ 1^ 0r (k) (k') (32) 

and the selfconsistency condition 

A(k) = -gr0r(k)i£#(k') {c^^j^cyj^) (33) 
k 

We adopt the Nambu formalism to diagonalise the mean-field Hamiltonian in eq. (22) which 
allows us to reduce it to single-particle form: 



AFBZ 

H= ^]|j{me(k)t3+hi+mA(k)fi}^,j (34) 

k 



where we explicitly exploited the fact that the (pseudo) spin is conserved. The Nambu spinors 
have 32 components and are defined as 

^k= (%T''^k+QiT'---''^-ki''^-k-Q,i''") ^^^^ 

They account for the coherent superposition of particles and holes which is the characteristic 
feature of the superconducting state. Here 1, fi and fa denote the unit matrix and the Pauli ma- 
trices in particle-hole space. The 16 wave vectors Qo = 0, Qi, • • • , Q15 are the reciprocal lattice 
vectors appearing in the antiferromagnetic phase. The set includes the eight propagation vectors 
of the SDW and their harmonics. The k-summation is restricted to the reduced Brillouin zone 
(AFBZ) of the antiferromagnetic state defined by the SDW. The structure of the Hamiltonian 
matrix in particle-hole space is 

me(k)+h mA \ 
' mA -me(k)+h J 
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where the 16 x 16-diagonal matrix contains the quasiparticle energies of the paramagnetic nor- 
mal phase 

™e(k)Q,.Q. = 5Q,.Q,8(k + Q0 (37) 

The magnetic Umklapp scattering associated with the modulated spin density is accounted for 
by the matrix 

h=-/jo(r)m (38) 

whose temperature-dependent amplitude ho {T) has to be determined selfconsistently while m 
is a given matrix. The diagonal matrix 

mA(k)Q,.Q. = Ao(r)m^>(k)Q,.Q. = 5Q,Q.Ao(r)0r(k + Q,) (39) 

contains the superconducting order parameters with the given k-dependent function 0r (k) and 
a temperature-dependent amplitude Aq{T). The self consistency equations in eq. (24) can be 
formulated in terms of the off-diagonal elements of the 32 x 32-matrix Green's function 

G{ien,k)^{ie„l-H{k))~^ (40) 
according to 

jj Ec AFBZ 1 

ho{T) = -r£ -Tr[miG(ie„,k)] 

i£n k 

Ec AFBZ 1 

HT) = -f^r£ £ -Tr[m«I>(k')fif3G(/e„,k')f3] (41) 

^ i£n k' 

Here /e„ = nT {In + 1) denote the T-dependent Matsubara frequencies and is the energy cut- 
off required in weak-coupling theory. The coupling constants U and as well as the cut-off 

Cc are eliminated in the usual way in favor of the decoupled transition temperatures 7^^'^'' and 

Tc'^\ The temperature dependence of SDW (A-phase) and SC order parameters is shown in the 
right panel of fig. 13. On entering the SC regime the A-phase is expelled in a finite temperature 
interval. This is nicely confirmed by recent neutron diffraction results [89] shown in the left 
panel of the figure. 



3.2 THE (ALMOST) IDEAL NFL COMPOUND CeNi2Ge2 

This is one of the few Ce-HF compounds which neither order magnetically nor become su- 
perconducting, except for incipient SC at very low temperatures in very clean samples. In- 
stead it exhibits clear cut NFL anomalies over more than two temperature decades. Together 
with the recently discovered YbRh2Si2, CeNi2Ge2 is one of the rare cases where one is ac- 
cidentally very close to a QCP at ambient pressure. This is supported further by the x (Pd 
-concentration) dependence of the Neel temperature in the alloy series Ce(Nii_;cPdx)2Ge2 [90]. 
It is found that Ta/(x) extrapolates to zero at at a very small positive concentration x^. Thus 
CeNi2Ge2 is (in contrast to YbRh2Si2) slightly on the nonmagnetic side of the QCP. After 
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Figure 13: Left panel: Integrated neutron diffraction intensity of the A-phase SDW satellite 
corresponding to Q ~ (0.215,0.215,0.530) as a function of temperature. It is proportional 
to the square of the A-phase order parameter in CeCu2Si2 [89]. Right panel: Variation with 
temperature of the A-phase (SDW) and superconducting order parameters with Fs-symmetry, 

^r-i (k) ~ coskxtt — cos kyU calculated for T^^^ =O.SK> Tc^^ = OJK [82, 83]. For this choice of 
parameters, the magnetic A-phase initially suppresses the superconducting state which forms at 
Tc < T^^^. The two ordering phenomena coexist in a narrow temperature range below Tg. The 
superconducting correlations, however, suppress the long-range magnetic order which disap- 
pears continuously at 7^ < < 7^^*^^ . 

some initial doubts [75] it is now accepted [91] that the 3D-SDW QCP scenario analysed in de- 
tail in [49] is realised in CeNi2Ge2. This predicts a scaling /(T) = 70 - cT^ which is well fulfilled 
in CeNi2Ge2 (fig. 14) down to temperatures around 0.3 K. For even lower temperatures sample 
dependent upturns in /(T) of unknown origin are present. Likewise the expected NFL resistivity 
behaviour p(T) = po -1- j8T^ with e > 1.5 depending on the sample quality (po) is found. The 
very closeness of the AF QCP implies an extreme sensitivity of the scaling to an external field. 
With increasing field strength one crosses rapidly into a Fermi liquid regime as witnessed by 
the appearance of a widening plateau in 7(T) (fig. 14) and a LFL type resistivity p(T) = po + 
A(B)T^. The boundary between the T*^ and T^ scaling marks the transition from the NFL to 
the LFL regime (fig. 14). As this boundary is approached from B > the A(B)-coefficient 
diverges indicating a singular energy dependence of the scattering rate. Thus as the QCP is ap- 
proached by varying B and T the quasiparticle mass and scattering rates in CeNi2Ge2 grow in 
agreement with the predictions of the 3D SDW scenario for quantum criticality. More recently 
the NFL behaviour in CeNi2Ge2 has has also been found in a low temperature divergence of the 
Griineisen ratio obtained from thermal expansion as discussed in sect. 2.3. 

3.3 SUPERCONDUCTIVITY AND QUANTUM CRITICALITY UNDER PRESSURE IN 
CePd2Si2 ANDCeRh2Si2 

At ambient pressure CeCu2Si2 and CeCu2(Sii_;iGe_^)2 together with CeCoIns and Celrlns are 
still the only Ce-based HE superconductors, aside from incipient superconductivity around 100 
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Figure 14: Left panel: Field dependence of 7 = C/T. Dotted line (B = 0) indicates y(T) = 70- 
Inset shows 7(B) in FL (plateau) regime [91]. Right panel: (a) B-T phase boundary 
separating the FL/NFL regimes; (b) Field dependence of the A-coefficient, solid line indicates 
A~B-°-6 [91]. 



mK in CeNi2Ge2 in exceptionally clean samples and the recently discovered non - centrosym- 
metric SC in CePtsSi [64]. In the former case this may be due to the accidental closeness of 
the A-SDW phase QCP at ambient pressure, witnessed by the phase diagram in fig. 5 and the 
associated NFL anomalies. To find more examples of Ce-based HF superconductors it is there- 
fore an obvious strategy to look for magnetically ordered Ce-HF compounds and drive them 
to a QCP with applied pressure, hoping that a SC 'dome', however small, might appear. This 
has been succesful for CeCu2Ge2 [92] and more recently for CeRh2Si2 [93] and also Celns and 
CePd2Si2 [94]. Note however that strictly speaking CeRh2Si2 is not a HF system due to its 
small mass enhancement visible from the low 7- value at ambient pressure. 

The magnetism of these compounds (table 2) and of the Pd-alloy series Ce(Rhi_;(^Pd;(^)2Si2 at 
ambient pressure has been subject to many investigations [94-96]. While there is unanimous 
agreement on the localised moment nature of CePd2Si2 , the interpretation of magnetism in 
CeRh2Si2 is controversial, partly favoring the itinerant picture from resistivity and specific heat 
measurements [97] and partly a localised picture from dHvA experiments and their interpre- 
tation and comparison with band structure calculations for the reference compound LaRh2Si2 
[98]. 

The commensurate AF structure of CePd2Si2 consists of FM (110) sheets stacked alternat- 
ingly with a wave vector Q = (^,^,0) and moments oriented along [110]. The local moment 
nature of CePd2Si2 can directly be infered from susceptibility measurements and INS experi- 
ments [99]. The former give a tetragonal CEF level scheme of the localised Ce^+ 4f^ electrons 
as a sequence of three Kramers doublets r^j\o), re(l9 meV) and ^(24 meV). Well developed 
AF spin waves with an uniaxial anisotropy gap A = 0.83 meV and a dispersion of 2 meV were 
observed which can be explained within the AF Heisenberg model for local moments. This is 
further supported by the magnetic phase diagram of the Ce(Rhi_;cPdjc)2Si2 alloy series shown 
in fig. 15 on the Pd-rich side: The monotonous decrease of Tyy and simultaneous strong increase 
of the characteristic HF temperature T* with Rh-doping (1-x) corresponds to the QCP of a Do- 
niach phase diagram for a Kondo-lattice type local moment system. There the destruction of AF 
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7(p=0) [mJ/molK^] 






T^(P«,) [K] 


Pm [GPa] 


CePd2Si2 


250 


10 


0.62 


0.4 


2.71 


CeRh2Si2 


23 


36 


1.42 


0.4 


1.05 



Table 2: Material parameters and magnetic properties of the sister compounds CeRh2Si2 and 
CePd2Si2. 




Figure 15: Magnetic x-T phase diagram 
of Ce(Rhi_;(;Pd;c)2Si2 with Neel temperature 
TnOO (circles) and Kondo temperature T*(x) 
(diamonds). Rh-rich side: Itinerant SDW 
regime. Pd-rich side: Local moment regime 
with Doniach type phase diagram [97]. 
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order is due to the compensation of magnetic moments in the Fermi sea of conduction electrons. 
The magnetism in Ce(Rhi_;cPd;c)2Si2 is less well understood. The origin of the exceptionally 
high Tivi = 36 K is a longstanding mystery which defies understanding in the local moment 
system if deGennes scaling of Ta^i with respect to GdRh2Si2 is applied. Like CePd2Si2 it or- 
ders with Q = (5,5,0), first with a single q-structure at Tjvi which changes to a 4q-structure 
below = 24 K (fig. 16). Also the (partly) plateau-like behaviour of Tiv(x) and T*(x) on the 
Rh-rich side of fig. 15 speaks against the local moment Kondo lattice picture and therefore an 
itinerant SDW origin of AF in Ce(Rhi_;cPd;c)2Si2 has been proposed in [97] . On the other hand 
the moment size (1.5/iB, m || c) can be explained within a CEF-split localised 4f^ model for 
Ce which is supported by the entropy gain AS(TAf) ~ Rln2 indicating a CEF Kramers doublet. 
From dHvA experiments [98] it was found that some Fermi surface branches are well explained 
by LDA results for LaRh2Si2 which advocates that 4f- electrons in CeRh2Si2 do not contribute 
to the Fermi surface. Alloying leads to an increase of the mean free path and may easily destroy 
unconventional superconducting states. Therefore, unlike in the CeCu2(Sii_xGe_^)2 system, no 
superconductivity has been found around the magnetic QCPs of the above T-x magnetic phase 
diagrams. However the magnetic QCP may also be approached by applying hydrostatic pres- 
sure to the pure compounds. Detailed studies of CePd2Si2 [95] have revealed the QCP to lie 
at pc = 2.86 GPa. The almost linear scaling Ta/(p) ~ (Pc-p) observed above 1.5 GPa is not 
in agreement with 3D spin fluctuation theory which predicts a pressure exponent 2/3. The re- 
sistivity at Pc exhibits a NFL-type dependence of p(T) = po+BT" with an exponent oc ~ 1.2. 
For samples with sufficiently small po an approximately symmetric superconducting dome was 
indeed identified around pc as shown in fig. 16. From upper critical field studies the coherence 
length was estimated as £,0 ~ 150A. 

Likewise hydrostatic pressure application destroys AF in CeRh2Si2 [100] at p^ = 1.05 GPa. 
The transition seems to be of first order. Nevertheless in good polycrystalline samples it was 
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Figure 16: Magnetic and SC phase boundaries under hydrostatic pressure in CePd2Si2 (left 
panel) [95] and CeRh2Si2 (right panel) [100]; note the enlarged scales of T^. In CeRh2Si2 the 
AF structure changes at T^2- 



found that SC appears around p^, albeit in a very small pressure interval from 0.97 to 1.20 
GPa, much smaller than previously thought (fig. 16). Again the optimal Tc ~ 0.4 K. From the 
Hc2 curves [98] at pc the coherence lengths — 240A and — 310A are obtained. The 
anisotropy of H^2^(0) may be caused by the anisotropic paramagnetic limiting effect {xJXa 
= 4), its observation also suggests a SC singlet pairing state. The dHvA results [100] show a 
discontinuous change of FS topology at p^, also the effective masses increase by almost an order 
of magnitude to m*/m ~ 30 above p^. Therefore the transition across the QCP was interpreted as 
change from the low pressure AF local moment phase with light electrons to the paramagnetic 
heavy electron state above pc [100]. 

Although the superficial appearance of quantum critical phase diagrams in fig. 16 is quite 
similar for both compounds, there is an important difference: quite opposite to CePd2Si2 , 
CeRh2Si2 does not show any signature of NFL behaviour around pc, y{T) has no significant 
T-depencence below 10 K and its value even increases up to 80 mJ/molK^ at pc due the sup- 
pression of AF order [101]. This means there is no direct connection between NFL behaviour 
in the QCP regime and the appearance of a SC dome. This conclusion can also be drawn from 
the investigation of Cel55 compounds (sect. 3.4). 

Indeed from the theoretical point of view there is no reason to expect such a connection as 
can be concluded from the discussion in sect. 4. At the QCP the magnetic correlation length 

diverges as described by eq. (9). Accordingly the spectral maximum at an energy ~ ^ 
of critical spin fluctuations responsible for NFL behaviour shifts to zero on approaching the 
QCP (eq.(45)). As explained in sect. 4, this increases the pair breaking effect due to low lying 
spin excitations and hence should lead even to a suppression of T^ close to the QCP. The soft 
spin fluctuations responsible for NFL behaviour in the Cel22 compounds close to the QCP are 
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therefore not essential for the SC pair formation. Therefore compounds like CeCu2Ge2 and 
CeRh2Si2 which are in the LFL state around Pc may form the same kind of SC state as NFL 
compounds like CePd2Si2. 

3.4 THE NEW Cel 15 AND Ce218 CLASS OF HE SUPERCONDUCTORS 

Only recently a new class of promising Ce-based HE superconductors with the general formula 
Ce„M,„In3„_|_2m (M = Co,Ir,Rh) has been discovered. Again the guiding principle of looking for 
SC close to magnetic QCP's has been succesful, which naturally means that for the members 
of the family the SC transition takes place out of a pronounced NFL-like normal state. These 
tetragonal compounds (space group P4/mmm) are composed of alternating n-fold Celns lay- 
ers and m-fold MIn2 layers derived from the parent compound Celns (n=oo). The n=l, m=l 
compounds (fig. 4) CeCoIns and Celrlns are ambient pressure superconductors and the AF 
CeRhlns becomes superconducting under pressure similar to the AF n=2, m=l family members 
Ce2CoIn8 and Ce2RhIn8. Although they partly have FS sheets with 2D character like slightly 
warped cylinders oriented along c, the physical properties, e.g. resistivity and upper critical 
field are not excessively anisotropic. In fact, an analysis of uniaxial pressure effects of both 
normal-state and SC properties based on high-resolution thermal expansion measurements re- 
veals that the SC T^ is strongly affected by at least two factors: The lattice anisotropy and the 
4f -conduction electron hybridisation which is most sensitive to c-axis lattice distortions [102]. 

Their electronic structure is well understood and the ambient pressure superconductors 
CeCoIns and Celrlns have conduction band states with strong 4f admixture, contrary to the 
ambient pressure AF compounds which have well localised 4f states witnessed by the similarity 
of their FS to the La parent compound. Large mass enhancements and specific heat 7 coeffi- 
cients are found but achieving the FL state may require application of a sufficiently large field. 
Maximum TcS ~ 2.5 K are achievable under pressure and already realised at ambient pressure 
for CeCoIns. This is higher than in other families of Ce- or U-HF superconductors, except for 
CeCu2Si2 under pressure, and a factor of ten larger than in the parent compound Celn^. A 
compilation of physical data is given in table 3, see also [103]. Certainly CeCoIns has turned 
out to be the most interesting system and is also most thoroughly investigated, partly due to its 
large ambient pressure Tc and large Hc2 as well as availability of excellent single crystals. Its 
unconventional superconducting state with d-wave symmetry is possibly the first which exhibits 
a change from a second to first order SC transition at Ht2(T) below a temperature Tq < T^ . Fur- 
thermore specific heat, NMR and ultrasonic attenuation suggest that CeCoIns is the first SC 
which exhibits the elusive Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase at high field and 
low temperature. 

3.4. 1 BASICS OF THE ELECTRONIC STRUCTURE 

There are extensive LDA calculations and dHvA investigations of the CeMIns and Ce2MIn8 
electronic structure and FS topology [100, 104—106]. For the reference compound LaRhlns ex- 
cellent agreement is found and the dHvA frequencies of a few sheets exhibit clear 2D cylindrical 
structure. The AF compound CeRhlns has localised 4f electrons which are CEF split into three 
Kramers doublets r^(0), rj(6.9 meV) and r6(23.6 meV) [107]. Therefore CeRhlnj should have 
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Table 3: Material parameters of Gel 15 and Ce218 HF superconductors. 



a similar FS as the La parent. Its AF order introduces many new sheets by folding into the AF 
BZ but the main FS sheets are indeed in good agreement. On the other hand the almost identical 
FS of CeCoIns and Celrlns are well explained by including the 4f states as itinerant electrons 
in the FS volume which makes them highly different from the above LaRhlns and CeRhlns FS. 
The selfconsistently calculated f-level occupations for CeCoIns and Celrlns are indeed close to 
one [105]. The enhanced y values are given in table 3 and are generally larger for the itinerant 
4f-electron compounds. The different size of mass enhancement for CeCoIns and Celrlns was 
attributed to the influence of CEF effects [105]. The mass anisotropies as obtained from upper 
critical field measurements [104] are considerable, m*/m* = 4.8 in Celrlns and m*/m* = 5.6 
in CeCoIns, reflecting the partly 2D FS topology. As in CeRhlns the AF bilayer compound 
Ce2RhIn8 has largely localised 4f states and therefore again the FS sheets are very similar to the 
La2RhIn8 parent compounds [106], three similar sheets with 2D cylindrical appearance have 
been found in dHvA experiments and LDA calculations. 

3.4.2 AF QUANTUM CRITICAL POINTS AND SUPERCONDUCTIVITY 

When magnetic order changes to a different symmetry or vanishes as a function of an exter- 
nal or internal parameter this may be interpreted as a quantum phase transition provided the 
respective ground state energies are largely determined by the contribution from quantum fluc- 
tuations. If the transition is of second order the critical parameter value then defines a proper 
magnetic quantum critical point (QCP). In its vicinity NFL behaviour and a superconducting 
phase transition may be induced by spin fluctuations of the type discussed in sect. 4. The Gel 15 
and Ce218 compounds give considerable support for this scenario. 

In the T-p phase diagrams of Gel 15 and Ce218 and also Gel 22 compounds which are AF 
at ambient pressure and SG under pressure two typical idealised situations may occur [108]: 
(1) When Tyv(p=0) ^ Tt(pt) a symmetric SG dome around p^ appears and FL behaviour is 
observed above T^ and for p > p^, AF order is suppressed inside the SG region. This situa- 
tion is realised for CeRh2Si2 and CePd2Si2 as seen in fig. 16. (2) When T;v(p=0) > Te(pe) 
the SG region is commonly asymmetric around the QCP and AF order partly coexists with SC. 
Around pc and above T^ an extended region of NFL behaviour is observed. This holds for 
GeGu2(Sii_xGe;i;)2, GeRhlns and GeRhi_;cG0;fIn5 (fig. 18). This picture is, however, oversim- 
plified as is obvious from the case of Gelns where the conditions (1) are fulfilled but nevertheless 
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coexistence SC/SDW has been identified in NMR experiments [79]. In fact these results suggest 
that instead of having a QCP within the SC dome a first order transition with AF/PM phase 
separation occurs. It remains to be seen whether this also holds true for CePd2Si2. Especially 



Figure 17: Left panel: Coexistence of AF order and SC in CeRhi-^CoxIns ('negative pressure' 
side) [109]. Right panel: Hydrostatic pressure dependence of in CeCoIns. The critical 
pressure where = is not known [100]. 

the quantum criticality of CeCoIns and its connection to superconductivity has been clarified 
by external field and pressure application and by internal pressure through Rh substitution of 
Co. At ambient pressure CeCoIns is nonmagnetic and the 'large' is almost independent in 
hydrostatic pressure up to 3 GPa and no magnetic order appears. Apparently doping with Rh, 
which has a larger ionic radius compared to Co, exerts a 'negative pressure' which for 1-x > 0.2 
Rh doping achieves AF order (fig. 17) while SC persists, albeit with slightly reduced Tf. Thus a 
considerable regime of AF/SC coexistence appears. Taken together Rh-doping and hydrostatic 
pressure do not present an ideal dome-shaped SC region, it is rather reminiscent of slightly 
doped CeCu2(Sii_xGex)2 under pressure. The combined x,p-T phase diagram of the Rh-doped 
compound CeRhi ;(^Ir;cIn5 [110] and Celrlns under hydrostatic pressure [111] is quite similar 
to fig. 17. CeRhi_;cIr;t:In5 also develops AF order that coexists with SC between 0.35< x <0.5 
[112]. The Tc(p) curve of Celrlns starts at Tc(p=0) = 0.4 K and reaches Tdpm) = 1.05 K at p^ 
= 2.2 GPa and finally drops to zero at pc = 6 GPa. This critical pressure for the CeCoIns sister 
compound is not yet known. 

Other family members which are AF at ambient pressure but not SC naturally exhibit better 
realisations of the canonical AF/SC quantum critical phase diagrams under hydrostatic pressure. 
This is already true for the AF parent compound Celn^ (T;v = lOK) of the Cell5 family [100] 
with an optimum T^ = 0.2 K at 2.5 GPa and especially for the localised 4f compound CeRhlns (x 
= in fig. 17) as is shown in fig. 18. The related AF bilayer compound Ce2RhIn8 is another 
exciting example of AF/SC quantum critical phase diagram (fig. 18). The remarkably high 
optimum T^. rivals that of ambient pressure CeCoIns and is comparable to its own T^(p=0) = 
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P^^fCOHj^j ^^^^^^ 

Figure 18: Quantum critical p-T phase diagrams for CeRhlns [100] and Ce2RhIn8 [113]. When 
^T'c in Ce2RhIn8, AF and SC coexist over considerable pressure range, contrary to the 
situation in CeRhlns where Tyy 3> T^. 



3.4.3 NFL ANOMALIES IN FIELD AND PRESSURE INDUCED QCP'S 

The common signature of a NFL state are anomalous temperature and field scaling exponents 
for thermodynamic quantities like specific heat, susceptibility, Griineisen parameters (thermal 
expansion) and transport quantities like resistivity or Hall coefficient. Their origin and appear- 
ance has been briefly discussed in sect. 2.3. 

Pronounced NFL anomalies have been identified in the ambient pressure SC CeCoIns and 
Celrlns, pointing to the closeness of the QCP [114]. This is well illustrated by p(T) data for 
CeCoIns at ambient pressure and zero field. They exhibit linear temperature behaviour above 
Tc and the logarithmic 7(T) variation between 0.4 K and 8 K for fields larger than H = H^^^^ = 5 
T but smaller than H = 8T (fig. 19) [62]. From systematic p{T,H) investigations of CeCoIns a 
B-T phase diagram may be mapped out showing the development of a field-tuned QCP and 
associated with it two different NFL regimes. They are characterised by p(r) = po+BT" where 
n=l or n=| and a low temperature crossover at T* to the LFL regime [1 15] is observed. The n=l 
regime is similar to the underdoped regime in high T^. superonductors and speculations about 
the existence of a spin pseudo-gap in CeCoIns and CeRhlns have been made [116]. This is 
supported by normal state Hall angle and magnetoresistance (MR) measurements [117] which 
show that in the n=l NFL regime the Hall angle &h ~ whereas Kohlers rule A^xiH) ~ 
which holds for LFL systems is strongly violated and a different scaling law is obeyed. The 
similarties of transport coefficients in CeCoIns to high T^ materials suggests a common origin 
of NFL behaviour in the critical 2D AF spin fluctuations in both type of compounds. This idea 
is also supported by an observed giant Nemst effect [1 18] as in the cuprates. Clear signature of 
NFL behaviour has also been seen in thermal expansion measurements [119] and the associated 
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Griineisen ratio. A compilation of various scaling laws in the NFL regime of the Cell5 and 
Ce218 compounds may be found in [103, 114]. A kind of inverted behaviour compared to 
CeCoIns shown in fig. 19 is seen in Cealring, because in this nonmagnetic HF compound LFL 
behaviour at zero field changes to NFL behaviour at a field of 13 T, which is interpreted as the 
appearance of a field induced magnetic QCP [120]. 

T(K) 

0.1 1 



CeCoInj 
*' J±H//c 

n=1 / n=2/3 
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SC ^ >^ ' n=2(LFU - 

5 10 15 
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Figure 19: Left panel: Field dependence of 7 = C/T for H || [001]. Solid lines are fits by 
Moriya's SCR spin fluctuation theory [62]. Right panel: B-T phase diagram of CeCoIns show- 
ing the evolution of two distinct NFL regimes (characterised by different exponents n of Ap(T) 
~ T«) and the high field LFL phase [1 15]. 




3.4.4 THE SUPERCONDUCTING GAP FUNCTION IN CeCoIns 

Heavy Fermion systems are unconventional superconductors regarding their pairing mechanism 
which is of magnetic origin in most cases and also with respect to non-s-wave symmetry of 
the order parameter, generally belonging to a non-trivial representation of the total symmetry 
group. Often, but not necessarily, this is connected with the appearance of nodes or zeroes 
in the gap function A(k). CeCoIns is the prime and best studied example of unconventional 
SC in the new family. It is a strong coupling SC with very large BCS ratio 2A/kTt = 8.86 at 
ambient pressure [121] and the gap function is believed to exhibit d-wave symmetry. This has 
already been suspected from the T^ -power law behaviour of NMR and NQR relaxation rate 
1/Ti and the absence of a Hebel-Slichter peak [122]. The same result was found for Celrlns. 
In addition CeCoIns shows a ^^^In- Knight shift reduction K ~ T for field along a and c, 
indicating an isotropic reduction of spin susceptibility in the SC state. The T- dependence of 
1/Ti and K consistently suggest a singlet pair state with line nodes, e.g. a d-wave state for 
CeCoIns. In subsequent field-angle resolved magnetothermal conductivity [60] and specific 
heat experiments [61] the position of line nodes, in k-space have been investigated to identify 
which irreducible D4/, representation among the five d-wave states is realised. In this method a 
peculiar feature of nodal superconductors in the vortex phase is used as has been first discussed 
by Volovik [123]. Because of A(k) = along certain directions in k-space, quasiparticles can 
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Figure 20: Specific heat y-coefficient for 
CeCoIns as function of field angle <j> with 
respect to [100] (H ± c-axis) in the ab 
plane. Maxima at [110] correspond to 
antinodal direction for A(k), suggesting a 
dxy-gap function. The oscillation ampli- 
tude decreases with increasing temperature 
[61]. 
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tunnel into the inter- vortex region where they acquire a Doppler shift energy Vs-k due to the 
superfluid velocity field Vi(r). This leads to a finite, field induced residual quasiparticle DOS. 
For a gap with line nodes one has (E <^ A): 

^^^|l^S^i((|£-v.(r).k|)) (42) 

where a double average over the Fermi surface and the vortex coordinate has to be taken. This 
contributes to the specific heat and most importantly supports a heat current ± to the vortices 
which is not possible in s-wave superconductors for T <^ T^. For H <^ Hc2 the residual DOS 
and hence 7(H) has a ^/H dependence characteristic of nodal superconductors [124]. In addition 
Ni(E,H) depends on the field direction with respect to crystal axes and the position of node lines. 
As a function of polar field angles 6 and (j) the residual DOS and hence the specific heat and 
thermal conductivity will exhibit angular oscillations whose type and phase allows one to draw 
direct conclusions on the k-space positions of nodes in A(k). The Doppler shift phenomenon 
has now been exploited succesfuUy to clarify the nodal gap structure in many unconventional 
superconductors [60, 125-128] including CeCoIns; see also [129] . There both Kxx(9,^) [60] 
and C(6,^) [61] exhibit fourfold oscillations in the azimuthal angle ^ which means that the gap 
function A(k) has four line nodes parallel to the c-axis on the main cylindrical FS of CeCoIns. 
The maximum in the oscillation will be achieved when each of them gives a Doppler shift con- 
tribution in eq. (42), i.e. when H points along the anti-nodal direction, the minimum occurs 
in the nodal direction. Surprisingly the two measurements presently disagree on the position of 
the line nodes: Thermal conductivity results suggest line nodes along c at [1 10] and equivalent 
directions. This means a d^2_y2 (Bi^) symmetry of the SC gap in CeCoIns as for high T^ com- 
pounds. On the other hand in the specific heat results of fig. 20 the maxima in the 0- oscillations 
are shifted by 45°, therefore node lines along c should be situated at the tetragonal [100] and 
equivalent positions meaning a dxy (B2g)-type symmetry of the SC gap. This discrepancy is not 
fully resolved so far. 
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3.4.5 VORTEX STATE AND FFLO PHASE IN CeCoIns 



Besides providing important clues on gap function symmetry the vortex state of CeCoIns has 
proved highly interesting, even unique, in its own right. It was known already for a long time 
that for low temperatures the SC transition at Hci should change from second to first order. 
This is due to a competition of orbital- and Pauli-pair breaking effects characterised by the 
Maki parameter a = \/2Hc2/Hjo where Hp = Ao/v^/is (for g = 2) is the Pauli limiting field 
. Conditions are favorable if the former is weak (large or a >1) compared to the latter 
and if spin-orbit scattering by impurities is negligible. Neglect of the orbital effect (a o°) 
leads to Tq = 0.55Tc for the appearance of the first order transition in an s-wave SC. Apparently 
these conditions have never been met in superconductors investigated so far. In CeCoIns, due 
to the large effective mass m* ~ 100 m and resulting extremely high H[y*^^ (0) = 1 1 .9 T, and due 
to excellent sample quality they are much better fulfilled. Indeed in this compound a change 
from second to first order superconducting transition for field H || [001] has been found at 
To — 0.3Tc (corresponding to a = 3.5) by specific heat [130, 131] and thermal expansion [102] 
measurements. The first order transition was also found for H || [110] from magnetisation 
[132] and specific heat [62] measurements below Tq = 0.5Tc. (fig 21). The same favorable 
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Figure 21: Left panel: Hc2 curve from specific heat results. At Tq — 1.1 K the SC transition 
changes from 1'^'^ to 2'"' order. Lower inset shows entropy gain. Upper inset shows specific heat 
jump at Tppj^Q line [62]. Right panel: enlarged FFLO corner determined from sound velocity 
measurements. Segmentation of vortices in FFLO phase with wave length 1/q is indicated [63]. 



circumstances may explain an even more exotic observation in CeCoIns for H ||[110]. At 
very small temperatures an additional phase appears in the vortex state (fig. 21). It is now 
considered as the first realisation of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase where 
the electrons form Cooper pairs with a finite momentum q. It appears when the Pauli effect 
largely dominates the orbital effect (oc > 1.8 [133]) and spin orbit scattering is very small. 
Apparently CeCoIns is the first SC compound where these conditions are fulfilled leading to 
the spatially modulated FFLO SC pair state with a wave number q ~ 2/igH/^v/7. In this state the 
order parameter should have planar nodes perpendicular to the field leading to a segmentation 
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of vortex lines along the vortex direction. A definite proof for this conjecture by small angle 
neutron scattering is still lacking, but indirect evidence for the segmentation has been found by 
sound velocity measurements [63]. The segmentation of vortices on crossing the Tpp^Q line 
changes their pinning properties which in turn leads to a Lorentz force anomaly for sound modes 
with displacement vector u ± H. The estimated segmentation length is remarkably close to the 
calculated modulation length 1/q = 60A. Also this method allows an independent determination 
of specific heat results and assures that the new phase is not due to field induced magnetic order 
since the latter should not have any effect on the vortex pinning. The FFLO comer is clearly 
visible for H in the ab- plane but hard to identify for H along [001] with the smaller Hc2 in 
agreement with theoretical predictions. Evidence for the FFLO state for H ||[110] has recently 
also been found from ^^^In-NMR measurements [134] which give a Tpp]^Q(H)-line similar to 
that from specific heat and sound velocity results and provide direct evidence for the spatial 
texture of the order parameter in the FFLO phase. 

3.5 HF SUPERCONDUCTIVITY WITHOUT INVERSION SYMMETRY IN CeRsSi 

All previously discussed Ce- intermetallic compounds and in fact all HF superconductors known 
possess centrosymmetric space groups which contain an inversion center. This leads one to a 
natural classification of SC pairs into even and odd parity or singlet and triplet states. For the 
singlet channnel only time reversal invariance is required in order to assure the necessary de- 
generacy of paired electron states with opposite momenta and spins. For the formation of triplet 
pairs however, electron states with opposite momenta but equal spins must also be degenerate 
which requires inversion symmetry as an additional necessary condition [135]. The lack of an 
inversion center may also strongly affect the magnetic properties of a superconductor as shown 
in [65]. 

The discovery of the first HF superconductor CePtsSi [64] without inversion symmetry is 
therefore of considerable fundamental interest. CePtsSi belongs to the tetragonal space group 
P4/mm (fig. 22) where the removal of the inversion center leads to the Ce-point group C^^ 
meaning the loss of the basal plane as mirror plane (z— >-z). The structure can be derived from 
the hypothetical cubic CePts structure (the same as for Celns) by tetragonal distortion with 
c/a=1.336 and by filling the voids with Si. The Ce^"'"4f states of CePtsSi are well localised 
as evident from sharp CEF excitations found in INS experiments [137] where the C4v level 
scheme was determined as T^^^O), r6(l meV) and ^24 meV). At T^ = 2.2 K AF magnetic 
order appears where FM planes are staggered along c with Q = (0,0,^). Below T^ = 0.75 
K superconductivity sets in (fig. 23) and coexists microscopically with AF order. Only the 
second time after A-type CeCu2Si2 (fig. 5) this has now been observed for a stoichiometric 
Ce-compound at ambient pressure. The ordered moments /i = 0.16jU5 however are much less 
than the value of the localised F^^^ ground state doublet which suggests a large Kondo screening, 
and indeed there is an associated mass enhancement of itinerant quasiparticles as is evident from 
7„(CePt3Si) = 390 mJ/mol as compared to y„ = 9 mJ/mol for LaPtaSi. A considerable 
residual 7^ = 180 mJ/mol and a AC/7„Tc = 0.25 much smaller than the BCS value 1.43 has 
been observed. This raises questions on the quality of the polycrystalline samples although the 
estimate from Wci indicates that one is not in the dirty limit. Later resistivity measurements on 
single crystals [138] in a field indeed have identified a coherence length ^ ~ 100 A 1 (mean 
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Figure 22: Left panel: Crystal structure of CePtaSi belongs to the tetragonal space group P4mm, 
the conventional unit cell (a = 4.072 A, c = 5.442 A) is indicated [64]. It lacks the reflection 
symmetry z -z. Right panel: Phase diagram of CePtsSii-xGe;^ from ambient pressure (mainly 
Tjv) and hydrostatic pressure results (mainly Tf). Here 1 % volume reduction corresponds to 1 .5 
GPa [136]. 



free path). There the almost isotropic upper critical fields has been determined as H^jCO) = 3.2 
T and H^!2(^) ~ ^-^ ^ -^^^ which is considerably larger than the Pauli-Clogston limiting 
field estimated by Hp ~A/v^/iB ~ 1.2 T for g = 2. The pressure depencence of T^ which 
vanishes around p^ = 1.5 GPa has still some uncertainty due to the broad appearance of the 
resistive transition [138]. Nevertheless a combined p(V)-T phase diagram from Ge-doping and 
hydrostatic pressure results is shown in fig. 22. Whether the AF T/^ and the ordered moment 
vanish continuously and the SC dome extends on both sides of the magnetic phase is presently 
not known. 

The observation of Hc2(0) ^ Hp would rather advocate for triplet pairing, this confronts 
one with two major problems: i) the lack of inversion symmetry should disfavor triplet pairing 
according to the above arguments, ii) coexistence with fully developed AF order should be 
destructive for triplet pairs. Proposals to avoid this contradiction have been made in theoretical 
investigations [66, 139] which include an explicit inversion symmetry breaking but time reversal 
conserving term in the single particle part of the BCS Hamiltionian. It has the form of a Rashba- 
type antisymmetric spin-orbit coupling for conduction electrons Cks which is of the form (ma 
= Pauli matrices) 

//p = a gk ■ mass'ci^Cks' (43) 

k.ss' 

where gk = -g k is an antisymmetric function which may be obtained from the periodic potential 
of Bloch states and is conveniently normalised to one on the Fermi surface. For finite a this 
term lifts the twofold spin degeneracy of conduction bands (but preserves the k, s ^ -k, -s 
symmetry). It leads to i) modifications for both singlet and triplet gap equations, ii) a mixing 
term for singlet and triplet gaps which, however, is of the order a/Cf <^ 1 and therefore can be 
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Figure 23: Left panel: Field dependence of C/T in polycrystalline CePtsSi with SC transition 
indicated by vertical lines [64]. Right panel: Upper critical field in single crystals from p [138]. 



neglected. Likewise reduction in the singlet case is of order (a/ep)^ <^ 1. The triplet case 
is more involved and provides important clues on the effect of Hp. For a = all triplet states of 
a a given odd / and described by a d(k) vector gap function have the same T^.. Turning on the I- 
symmetry breaking (a >0) lifts this degeneracy, the expected Tc- reduction however does not 
affect all triplet states indiscriminantly, those triplet states with d(k)|| gk are well protected and 
have no T^ reduction as in the singlet case. From the €4^ point group one may expect gk ■ mcr 
= kxOy-kyOx and then the protected triplet state with unaffected T^ corresponds to the A2„ gap 
function d(k) = xky-ykx which has point nodes . Other triplet states are strongly suppressed 
for a:$>Tc. Thus in principle favorable triplet states are not destroyed by inversion symmetry 
breaking. 

On the other hand the loss of inversion symmetry described by eq. (43) weakens the Pauli- 
Clogston limiting argument against singlet pairing considerably. It is known that spin-orbit 
impurity scattering reduces the effect of paramagnetic limiting, the same is true for the anti- 
symmetric spin-orbit coupling of paired conduction electrons in eq.( 43). With increasing a the 
renormalised Hp increases above its starting value A/ A/2/iB for a = 0. For a/Tc S> 1 it becomes 
arbitrarily large. As a result the observation of a large Hf2(0) ^ A/ \/ljlB does not contradict 
singlet pairing, which in addition would be more compatible with the coexisting AF order. 

As a preliminary conclusion one may say that there is no compelling theoretical reason why 
triplet pairing for SC without inversion symmetry should be prohibited. Likewise the upper 
critical field observations do not contradict singlet pairing so that the question of gap symmetry 
in CePtsSi remains open. Recent NMR relaxation experiments [140] on Pt-sites also did not 
clarify this question but rather added new puzzles: The low temperature 1/Ti behaviour is 



40 



neither described by the usual T-^ power law generally observed in HF SC, nor by an exponential 
decay indicative of an isotropic gap. However the latter is more compatible with the observation 
of a small coherence peak below and a Balian-Werthamer type relaxation rate immediately 
below it. 



4 MICROSCOPIC SC PAIRING MECHANISM IN HF COMPOUNDS 

In this section we illustrate to some detail microscopic concepts and theories that were proposed 
to describe the superconducting pairing in HF systems or, more generally, in strongly correlated 
electron compounds including the cuprates and ruthenates. The long history of the subject of 
non-phononic superconductivity will not be recounted and we refer to the literature for this 
purpose [6, 141, 142]. 

Far below the Kondo temperature T* the ideal HF compound is in a LFL state with quasipar- 
ticles that have an effective mass m* much enhanced compared to the band mass m^,. However 
due to the presence of virtual high energy excitations a residual on-site quasiparticle repulsion 
appears. The low energy excitations may then be described by a spectrum of spin fluctuations 
characterised by a dynamic susceptibility matrix in;^(q, (o). In the vicinity of a magnetic in- 
stability, for example AF order close to a QCP, overdamped critical spin fluctuations called 
antiferromagnetic paramagnons appear whose spectral density is concentrated around the AF 
Q-vector and which has a typical energy Tv/^ < T* at the maximum of m;^(Q, co). The virtual 
exchange of these excitations mediates an effective interaction between quasiparticles which 
may lead to a superconducting instability. In contrast to conventional phonon-mediated super- 
conductors the pair states have to be of unconventional (non s-wave) type whose order parameter 
exhibits a sign change in k-space to avoid the strong on-site repulsion of quasiparticles. A brief 
symmetry classification of unconventional SC order parameters was given in sect. 3.1. 

Originally the spin fluctuation mechanism in the case of nearly FM systems with param- 
agnon exchange was proposed to explain superfluidity in ■'He [143, 144]. It was later extended 
to HF [145] and cuprate superconductors [146]. In the latter one has simple tight binding (TB) 
hole bands which allows one to carry out fully microscopic calculations based on the fluctua- 
tion exchange (FLEX) approximation [147, 148] for the Hubbard model. However the theory 
in its original single band form is not able to describe the complications of SC pairing in real 
HF materials, for the following reasons: i) Orbital (spin-orbit) and crystal field effects lead to 
complicated renormalised band structures with multisheeted FS, therefore the dynamic suscep- 
tibility mx (q, (d) cannot be reliably calculated, instead empirical forms have to be assumed 
which also incorporate the vicinity to a possible QCP frequently present in Ce-HF compounds, 
ii) In the U-compounds 5f- electron have a dual nature, partly localised and partly itinerant. The 
internal CEF excitations (magnetic excitons) of the former may dominate the effective interac- 
tions between itinerant quasiparticles, which is an alternative to the spin fluctuation mechanism 
that involves only the itinerant 5f electrons. 

We shall discuss both electronic pairing mechanisms in the following. While the spin fluctu- 
ation mechanism may be appropriate for the Kondo-like Ce HF superconductors, the magnetic 
exciton mechanism is a serious candidate for a number of U-compounds. 
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4. 1 SPIN FLUCTUATION MECHANISM FOR Ce-BASED HF COMPOUNDS 



If one is not too close to the QCP of magnetic order one may assume that the frequency depen- 
dence of spin fluctuations is negligible on the scale of the cutoff energy for SC pairing. Then in 
the effective interaction V(q,£0) between quasiparticles m;^(q, to) may be replaced by the non- 
retarded RPA expression mx{q). For the spin dependent part one obtains the spin rotationally 
invariant interaction 

^aj3;yp (q) = V(q)mc7ay • mcT^p with 



Here I is the residual on-site quasiparticle interaction and Xo{^) ^ their static noninteracting 
susceptibility. The interaction may be split into singlet (S=0) and triplet (S=l) channels accord- 
ing to Vo = -3V and Vi = V, the factor 3 is due to the fact that only longitudinal spin fluctuations 
contribute in the triplet channel whereas both longitudinal and transverse parts contribute in the 
singlet channel. In this approximation the SC pairing is completely determined by the static 
susceptibility m;^(q) which may in principle be obtained from diffuse neutron scattering. Its 
q- dependence together with the characterisics of the Fermi surface determines the dominant 
angular momentum component of Ysi^) that leads to the stable SC pair state (/ = even for 
S=0 and odd for S=l). Obviously the s-wave state is not stable because Vq"^ is always positive 
(repulsive), pair states with / > 0, however, may have attractive V/ < and become stable de- 
pending on the characteristics of x{*i)- The model has been analysed in some detail in [145] by 
assuming a spherical FS and a plausible behaviour for V(q) in a cubic lattice with a peak near 
the AF wave vector Q to simulate the AF spin fluctuations observed in HF compounds. This 
analysis predicts that the unconventional even singlet SC state (/ > 0), e.g. / = 2 d-wave state, 
should be stable because V/ < is attractive in this orbital channel. 

This result led to the general expectation that AF spin fluctuations favor singlet pairing. 
While it is true in the simple case of cuprates, and it may be true for some Ce-based HF com- 
pounds like CeCoIns (sect. 3.4), it turned out to be misleading for a number of U-compounds 
like UPt3 or UNi2Al3 and possibly UBei3 which exhibit triplet pairing as concluded e.g. from 
Knight shift results [149]. Therefore the above prediction of singlet pairing may be partly an 
artefact of the model simplifications and point to the inadequacy of the spin fluctuation model 
for U-HF superconductors. Indeed the model has obvious deficiencies: i) It completely ignores 
the orbital structure strongly influenced by spin orbit coupling and CEF effects, ii) It does not 
take into account the uniaxial crystal anisotropy present in most HF superconductors which 
possess either tetragonal or hexagonal space groups, iii) It ignores the spin space anisotropy of 
the pair potential which is assumed as rotationally invariant in eq.(44). In case this symmetry is 
reduced the triplet state will be less disfavored. 

From the low temperature 'power laws' in various thermodynamic and spectroscopic quan- 
tities it is also concluded that most HF SC have gap functions with nodes A(k) = along lines 
on the FS. However if one assumes the strong spin orbit coupling case for quasiparticles, i.e. the 
crystal double group as symmetry group, 'Blount's theorem' states that triplet pair states can 
have only point nodes in the SC gap [150]. Therefore the observed power laws which support 
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line nodes have also been used as an argument for singlet pairing. This argument is misleading 
however, since the spin orbit coupling is so large that it has to be included already in the orbital 
basis for band structure calculations within the jj- coupling scheme. The scale of SC pairing 
energies is much smaller than band energies and therefore it is the 'pseudo-spin' degrees of 
freedom connected with Kramers degeneracy of quasiparticle band states which are involved in 
the singlet and triplet pair wave functions. Therefore the pseudo-spins may be assumed to have 
an effectively weak (pseudo-)spin orbit coupling, and then no contradiction between the exis- 
tence of line nodes and SC pairing in triplet states for U-compounds appears. We note however 
that ironically, in UPdjAls, which is the one U compound where the spin fluctuation model is 
known to be inappropriate probably a singlet SC pair state is realised. 

In a series of papers the spin fluctuation mechanism has also been studied within the non- 
retarded strong-coupling Eliashberg approach [151-154] including the frequency dependence 
of the pairing interaction and assuming a simple tight binding band with n.n. and n.n.n hopping 
t and t' respectively for conduction electrons. The frequency dependence is especially important 
close to a magnetic QCP when the typical energy of spin fluctuations becomes comparable to 
InTc- Instead of starting from a microscopic model Hamiltionian, a phenomenological effective 
interaction kernel V(q,(0) was used as input for the Eliashberg equation which determines the 
quasiparticle self energy E(k,ift)n) and anomalous self energy ('gap function') A(k,i(On). For 
singlet (S=0) and triplet (S=l) pairing it is given in the AF case by 

^^(q, CO) = as^xifl, <») with ao = -1 and ai = ^ 

2 

^ = ^2 ^^2'^- 0) where 77 (4) = T,fq_ , q^q+ (45) 

q\ = (4-|-2am) ±2[cos^;c-|-cosq^-y-|- o^cosq^^] 

Here k and Kq denote the nearly critical and background inverse magnetic correlation lengths 
respectively (in units of the inverse lattice spacing a^). In Stoner-RPA theory = (1- 

lXo)~^ and Kq ~ Ikpa. This interaction is mediated by overdamped spin fluctuations since 
the pole of X{(i,0}) is at a purely imaginary frequency, i.e. the spin excitation does not have a 
real group velocity and therefore is not propagating. This is a natural assumption close to the 
QCP where the restoring force for spin excitations due to the molecular field breaks down. It is 
completely different from the phonon case and also the magnetic exciton mechanism discussed 
later which are both characterised by the exchange of propagating bosons which have a real 
frequency bigger than the line width. In eq.( 45) determines the crystal anisotropy of the 
spectrum, it is of 2D character when = and 3D character for = 1. It is still fully 
isotropic in (pseudo-) spin space as evident from the singlet and triplet classification. According 
to eq.( 45) the interaction in the AF case peaks at the zone boundaray vector Q with a spread 
given by k{k = inverse magnetic correlation length t,~^) and a frequency maximum at (Omax = 
77 (q) K^. The same type of expression holds for EM spin fluctuations where one has to replace 
q = q-. Using the effective pairing in eq. (45) the Eliashberg equation for the gap function 
A(k, ifflb) {(On = Matsubara frequencies) is 

T 

A(r) A(k, ia)„) = as - £ £ Z (k - q, i<On - i<Om) I G(q, i(0^) \ i(Om) (46) 
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where G(q,ift),„) is the quasiparticle Green's function renormalised by the normal self energy 
E(q,ift)m) that satisfies a similar equation. The solution of this set of self consistent equations 
determines Tc via the condition that the eigenvalue A(Tc) = 1 . The fixed parameters of the model 
are t,t' (hopping elements), Kq (inverse background correlation length), and T^/ (spin fluctuation 
energy scale). The critical temperature Tc(A,Jc) is studied as function of the dimensionless 
coupling constant A = l^Xo^^ the inverse magnetic correlation lenght K which determines 
the width of the AF spin fluctuation peak in k-space. In reality this parameter may be varied 
experimentally by applying pressure, i.e. tuning the distance to the magnetic QCP where K = 
0. In the spirit of Ginzburg-Landau theory, at T^ the gap function should belong to a single 
representation T of the spatial symmetry group whose k-dependence is characterised by the 
basis function ^•r(k). The stable SC state is the one which has the highest T^- among the 
representations T and pair spins S =0,1. The Tt(A,K') variation of the stable state may then be 
calculated. This has been done e.g. for the d^2_y2- state in the 2D case appropriate for cuprates, 
or, with suitable redefinition of bandwidth and T^/ also for CeCoIns. Without going to great 
detail the essential results with the above procedure can be summarised as follows [151-153]: 

• When Q = (FM) the p-wave state (S=l) is favored (has the largest T^) and for the AF 
Q = (^,i) (2D) or Q = (5,5,3) (3D) the d-wave state (S=0) is prefered. 

• In both cases T^ increases monotonously with coupling constant A for any K. For constant 
A it first increases with K and reaches a maximum at ~ 0.5 and then shows slow 
monotonic decrease with K. When A decreases the optimum T^. position moves towards 
JC = 0, i.e. to the AF QCP. 

• Generally the achievable Tc's under equal conditions are larger in 2D than in 3D, and they 
are much larger for AF spin fluctuations leading to even pair states as compared to FM 
spin fluctuations associated with the odd SC states. The latter also depend more strongly 
on parameters, e.g. on JC. 

This result is somewhat counter-intuitive: In the FM case the effective interaction in real space 
is purely attractive while it is oscillating in the AF case. However as mentioned above, for triplet 
pairing only longitudinal spin fluctuations contribute, leading to the reduction factor ai=l/3 in 
eq. (45) whereas in the AF singlet pairing case the additional transverse contributions lead to 
|(3o|=l- This advantage of singlet pairing depends on the isotropy of the effective interaction 
in spin space. In the opposite limit with maximal Ising type anisotropy in spin space, the 
disadvantage of triplet pairing is removed as will be discussed in sect. 5.2. 

4.2 Tc-DEPENDENCE ON SPECTRAL PROPERTIES OF SPIN FLUCTUATIONS 

In the above considerations the change of T^ under variation of global interaction parameters 
and characteristics is considered. This does not answer the question how the variation of the 
spectrum of spin fluctuations for a specific q* and (Oq influences T^, specifically one would like 
to know which parts of the spectrum are most favorable for achieving a large T^- for a given pair 
state. For the electron-phonon superconductor this has been first investigated in [155]. It was 
found that every part of the e-ph interaction function oc^F((0) contributes positively with 2kTc 
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Figure 24: Dependence of AT^- on spin fluc- 
tuation frequency (Ho (t = hopping energy) for 
various momenta q* in the 2D AF case with 
Q = (0.5,0.5) (in r.l.u.). Coupling constant 
^X^k = 10; = 2 = 0.5, K-2 = 12 and T,/ = 
0.67t. Only spin fluctuations at q* = (;r/a,;r/a) 
= Q or close to Q contribute positively to T^. 
above the crossover frequency (Ocr(q*) — 0.5t. 
Below C!)cr(q*) they are pair breaking. [154]. 



being the optimum phonon frequency while phonons with frequencies much larger or much 
smaller frequencies contribute little to T^. 

This analysis has been extended later [156] to the spin fluctuation model for HF systems, 
albeit using an unrealistic effective interaction that factorises with respect to frequency and mo- 
mentum. As a main result it was found that, for any momentum q, adding spectral weight below 
a crossover frequency cOcr reduces T^, i.e. low energy spin fluctuations act pair breaking. How- 
ever a subsequent analysis by McHale and Monthoux [154] using the more realistic effective 
interaction in eq. (45) has modified this picture considerably. By adding an infinitesimal amount 
of spectral weight to ;^(q, to) at fixed wave vector q* and frequency COq and -(Oo, i-C- adding two 
delta-functions, the change in the critical temperature Ar^ (q*, (Oq) was computed numerically. 
A representative result in 2D for three different q* is shown in fig. 24 for the nearly AF spin 
fluctuation spectrum with 2D incipient ordering vector Q = (i,i). It clearly shows that the effect 
of the added intensity to the spectrum also depends strongly on q*, i.e. where it is added in the 
BZ, in addition to the (Oq dependence. For q* = (0,0) or (^,0) which are far away from the AF 
Q = (^,^) the added intensity reduces T^- for all frequencies cOq. Only when it is added at the 
incipient ordering wave vector Q or close to it one has ATc > for frequencies CO > (Ocrin*)- 
The crossover frequency cOcr from pair breaking to pair formation therefore strongly depends on 
the wave vector q* of the added intensity relative to Q. Sufficiently far away from their spectral 
maximum at Q spin fluctuations with all frequencies (Oq are pair breaking (A71(q*, COq) < 0). 
The pair forming fluctuations are constrained to the region around Q which means that the 
pairing occurs mostly through the short range magnetic fluctuations, except very close to the 
QCP (k — > 0). The analysis again shows that ATc{q*, (Oq) is much more robust in the AF case 
as compared to FM case (Q = 0). From the Arc(q*, (Oq) curves in the AF case it is also found 
that the optimum coq where ATc is largest lies above the spectral maximum of x{fl, 0)) for <C 
1 (nearly critical case) and vice versa for ^ 1 . Despite the complications the optimum (Oq is 
of the same order of magnitude given by 2;rTc as for the electron-phonon mechanism. 

These results show that the formation of the SC state is not determined by the very low 
energy critical spin fluctuations which are responsible for the NFL anomalies in compounds 
like CePd2Si2. It explains qualitatively why compounds like CeCu2Ge2 or CeRh2Si2 which 
exhibit LFL behaviour around the critical pressure for destruction of AF order may nevertheless 
be spin fluctuation mediated superconductors (see also sect. 3.3). 
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The previous analysis also suggests another reason why triplet pairing should be relatively 
scarce. In strongly correlated elecron systems not too close to a QCP dominant spin fluctuations 
involve mostly nearest neighbors, i.e. they have wave vectors Q at the zone boundary, these are 
however triplet-pair breaking for all frequencies in the spin fluctuation mechanism. Considering 
this fact and the observation that most U-based HF superconductors seem to have triplet pairs it 
seems doubtful whether the spin fluctuation model is appropriate for them and one has to think 
about other possibilities. 



4.3 MAGNETIC EXCITON MECHANISM FOR SUPERCONDUCTIVITY IN U- BASED 
HF COMPOUNDS 

The origin of mass enhancement in Ce- and U-based compounds is not the same as explained in 
sect. 2. This is due to the multiple 5f-shell occupation and strongly orbital dependent hybridi- 
sation in U compounds which, in some cases like UPd2Al3, UPt3 and possibly UBei3 leads to 
a partly itinerant and partly localised nature of 5f electrons. In a strong-coupling theory for the 
quasiparticles, mass enhancement and effective pair interactions are two sides of the same coin. 
Therefore we discuss now the origin of SC pairing within the dual model of 5f electrons. This 
idea has been proposed and explored to some extent in [157-159] for UPd2Al3 as discussed 
further in sect. 5.2. 

In the dual model the localised 5f^ electrons occupy specific total angular momentum or- 
bitals which are split by the CEF. For simplicity we consider only the two lowest CEF states 
which are assumed to be nonmagnetic singlets. They are then described by a pseudo-spin (S=^) 
variable where Sj; = correspond to singlet CEF states with energies ±|. They have an effec- 
tive inter- site coupling (superexchange) as well as a coupling to the itinerant 5f electrons of the 
conduction band e^. The latter couples only with the o^- conduction electron spin component 
due to the (Ising-type) anisotropy introduced by the singlet CEF states which have only transi- 
tion matrix elements for the local S^ pseudo-spin operator. This dual 5f-model is described by 
[159] 

^ = L ekcL^kcT + 5 £ - 7// £ S^Si - 1 £ Gi,Si,. (47) 

k(T (• ((7) i 

The localised part may be diagonalised separately. Due to the effective inter-site interaction 
J(q) which contains the superexchange J// and a RKKY contribution from the last term the 
local CEF excitations at an energy 5 evolve into propagating magnetic exciton modes with a 
dispersion given by 

coEiq) = 5[l-:^tanh|5] (48) 

If J(q) has its maximum Jg = J(Q) at an AF zone boundary vector Q the mode frequency be- 
comes soft at Q at the Neel temperature where induced AF appears; 

provided one has > 1 for the control parameter. In reality the softening will be arrested at 
a finite magnetic exciton gap due to the effect of higher lying CEF states which may shift T^v 
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slightly to higher temperatures. The magnetic excitons associated with the localised 5f system 
may be seen in INS as discussed for UPd2Al3 in sect. 5.2. For T ^ 5 they are bosonic modes 
like phonons and then a canonical transformation on eq. (47) leads to a more convenient dual 
model Hamiltonian [159] 



where c^^ and ttq create conduction electrons and magnetic excitons respectively and V^a(r), 
^ (r) are the associated field operators. It is indeed very similar to the electron-phonon Hamilto- 
nian, except that magnetic excitons do not couple to the electronic density but (in this simplified 
model) only to the z-component of the itinerant spin density. The resulting normal self energy 
of conduction electrons due to dressing by virtual magnetic excitons has been discussed before. 
Complementary virtual exchange processes involving magnetic excitons lead to an effective 
quasiparticle interaction given by (v„ = 2;rnT) [159] 



This effective pairing breaks spin rotational symmetry in a maximal (Ising type) manner. This 
is enforced by the anistropy of the CEF states and has important consequences for the classifica- 
tion of SC pair states, one now has to use the equal spin pairing (ESP) states \x) — |TT) and 
and opposite spin pairing (OS?) states, given by 



This aspect is completely different from the previous spin fluctuation theory which leads to a 

spin rotationally invariant pairing potential and the pair states are classified as (pseudo-) spin 
singlet (S = 0) and triplet (S = 1) states. As a consequence, because p ={x\^^^^\x) is +1 for 
ESP and -1 for OSP respectively, the coupling strength of the effective potential in eq. (51) is 
equal for OSP and ESP states, in contrast to the enhancement factor 3 of singlet vs. triplet pair 
interaction in the spin fluctuation model. Thus the present mechanism with its Ising type spin 
anisotropy in eq. (51) does not disfavor the odd parity states over even parity states, in fact they 
may be degenerate as will be discussed for a concrete case in UPd2Al3. 

The magnetic exciton mediated pairing mechanism is due to the interaction between itiner- 
ant 5f electrons caused by the exchange of excitations within the localised 5f-CEF states. On 
one hand it is similar to the electron-phonon mechanism because propagating modes are in- 
volved and not overdamped spin fluctuations, on the other hand it is repulsive in the s-wave 
channel and contains the conduction electron spin variables leading to the possibility of uncon- 
ventional or nodal pair states. 



5 U- BASED HF SUPERCONDUCTORS 

Most ambient pressure HF supercoductors are intermetallic U compounds like UBeo [160], 
UPt3 [161], URu2Si2 [162], UPd2Al3 [163] and UNi2Al3 [164]. This may partly be caused by 




kcr q 



(50) 




(51) 




(52) 
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the enhanced delocalisation of 5f as compared to 4f electrons. Due to the multiply occupied 5f 
shell in U which hosts between two and three electrons, the origin of heavy electron mass en- 
hancement in U compounds is in fact quite different as compared to Ce compounds with nearly 
singly occupied 4f orbitals. In the latter, heavy quasiparticles originate from coherent resonant 
scattering in the Kendo lattice [4]. For some U-HF systems however the dual, i.e. localised and 
itinerant nature of 5f electrons in different orbitals generates the mass enhancement [28, 29] 
(sect. 2) and may also dominate the SC pairing mechanism. The AF QCP scenario so prevalent 
in Ce-compounds is not of general importance in U-based HF superconductors with the possible 
exception of UBei3. One rather observes that SC phases in U-compounds are embedded in a 
region of reduced moment AF order. In the more recently discovered FM U-based superconduc- 
tors like UGe2 or URhGe however QCP of FM or even hidden order may indeed be associated 
with the observed unconventional SC phases. On the other hand these superconductors are not 
really HF materials as witnessed by their low 7- values and they will not be discussed here. 

In this review we shall not give a complete overview on U-based HF superconductors but 
shall focus on two examples, Ui_;cTh;cBei3 and UPd2Al3 where new results have been obtained 
recently. Some of the other classical U-compounds, notably the unique multicomponent super- 
conductor UPt3 are reviewed in [6, 8, 9] while URu2Si2 is reviewed in [6]. 



5.1 MULTIPHASE SUPERCONDUCTIVITY IN THE NON-FERMI LIQUID COMPOUND 

Ui_;cTh;,Bei3 

The cubic Be-cage compound UBeo was the second HF superconductor discovered [160, 165] 
but remains only partly understood. The SC transition takes place in a normal state with pro- 
nounced NFL behaviour, as can be seen from the C/T behaviour for fields above Hc2. The 
thorated UBei3 has a complicated x-T phase diagram with two distinct SC phases and possibly 
a coexisting magnetic phase. In fact it is surprising that SC survives up to about x = 0.06 since 
the potential scattering of normal Th impurities should lead to a strong pair breaking for an 
unconventional SC state and an ensuing a rapid reduction of Tc. Instead a pronounced non- 
monotonic behaviour of Tc(x) is found. The x-T phase boundaries are still partly a conjecture 
and may have to be revised according to recent flux creep experiments [166]. The most interest- 
ing issue is doubtlessly the origin of small magnetic moments in the low temperature SC phase 
(C) in the regime 0.02 < x < 0.045. They may either be caused by a non-unitary SC C-phase 
which breaks time reversal symmetry as suggested by juSR experiments [167] or they may be 
associated with the formation of a long range SDW order suggested by thermal expansion ex- 
periments [168-170]. Theoretically the former scenario has been explored within the context 
of a Ginzburg-Landau approach assuming that the Tt(x) of two distinct SC phases cross [5]. 
The SC/SDW scenario on the other hand has been theoretically proposed in [171]. It should be 
noted however that assumptions on the gap function symmetry remain purely speculative. Even 
in the pure UBei3 it is not known although the presence of nodes in the gap function and hence 
its unconventional nature is assured, e.g. from the existence of a zero-bias conductance peak 
in tunneling experiments [172] which provides evidence for a sign change of the gap function 
A(k). 
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5.1.1 ELECTRONIC STRUCTURE AND NFL-LIKE NORMAL STATE 



The 5f states in UBei3 show both signatures of localisation and itineracy. On one hand the 
5f- orbital energy is close to the Fermi level according to photoemission results. This disfavors 
the Kondo picture with a localised 5f^ configuration as origin of the extremely large mass en- 
hancement characterised by 7= 1100 mJ/mol corresponding to a quasiparticle band width 
of T* ~ 8 K. LDA calculations [173] with fully itinerant 5f electrons however lead to only 7, 
= 13 mJ/mol implying a large many body enhancement factor of m*/m^ = 85. The Fermi 
surface consists of two closed hole sheets with nesting properties and one closed electron sheet. 

On the other hand Schottky-type specific heat anomalies [174] suggest the presence of CEF 
excitations ~ 180 K out of a r6(5f^) Kramers doublet ground state. An alternative model [175] 
assuming a nonmagnetic r3(5f^) ground state can be ruled out because it predicts the wrong 
sign of the nonlinear suscetibility X^^\T) low temperatures [176]. A mixed valent model 
involving both nonmagnetic and magnetic Eg doublets has also lead to the conclusion that 
the ground state must be the magnetic Eg to obtain qualitative agreement with the observed 
;^*^-^)(T) behaviour [177]. In the pure Eg picture the Kondo effect would lead to a LFL state 
with enhanced m*. However a further temperature scale below T* given by T^ ~ 2 K has 
been found from a clear maximum structure in resistivity, specific heat and thermal expansion 
[169, 170]. This indicates that UBei3 has not yet reached a true LFL state when SC sets in at 
Tc ~ 0.9 K. Indeed for fields above Hci a pronounced NFL behaviour with logarithmic increase 
in y(r) and no saturation to the lowest temperatures has been observed [178, 179] in UBeis. 
In addition a universal field-independent p{T)/piK behaviour is observed for 7 > 1.2^ and 
AT <H < lOr in the normal state. One finds p{T) r^T at elevated and (p(r) - po) ~ T^/^ 
at lower T [179, 180]. The origin of NFL behaviour may reside in the vicinity of a QCP 
for SDW or possibly hidden order, this would be in accordance with the nesting properties of 
electronic bands in UBeis and with the evolution of TmOO which hits the maximum of Tc(x) 
in Ui_xTh;cBei3 as concluded from resistivity measurements [169, 170]. Another possibility is 
provided by the T^(B*) dependence, implying T^ — > for B* — > 57 (fig. 25), i.e. a field-induced 
QCP of short-range AF correlations [179]. 

5. 1 .2 THE SUPERCONDUCTING PHASE DIAGRAM OF Ui_;,Th;cBei3 

The symmetry of the SC state in pure UBeo is not known with any confidence, but its anisotropic 
nature is suggested by power law behaviour of specific heat, penetration depth and NMR relax- 
ation and by tunneling results mentioned before. As usual from the power law exponents no 
consenus has emerged whether the gap function is characterised by point or line nodes and 
therefore we do not discuss these results. A strong indication of nonconventional SC is ob- 
tained from a giant ultrasonic absorption anomaly immediately below Tc [181, 182] which was 
attributed to damping by domain walls formed by a multicomponent SC order parameter [5]. 
Also the B-T phase diagram (fig. 25) is quite anomalous, it exhibits a line of anomaly B*(T) 
deep in the SC regime which starts at Tl(x=0) = 0.7 K signifying the onset of short range AF 
correlations. As function of x they eventually develop into a long range SDW order discussed 
below. The SC phase diagram of thoriated Ui_;cTh;cBei3 (x < 0.1) has required a great effort 
to unravel and the task is still not finished. It displays a bewildering richnes of different phases 
and lines of anomalies which have been obtained mostly from specific heat, resistivity and ther- 
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Figure 25: Left panel: B-T phase diagram of UBei3 with additional anomaly line B*(T) ob- 
tained from specific heat (squares) and thermal expansion (triangles) experiments. Right panel: 
SC phase diagram of Ui_;(;Th;cBei3 . Full lines and symbols denote boundaries of different ther- 
modynamic phases, broken lines and open symbols denote lines of anomalies. Tl and T^-^ are 
anomaly lines obtained from the minimum in thermal expansion a(T). = is the line of 
anomalies from the maximum in C(T) and a(T). A,B and C are the distinct SC phases in the 
Tc-crossing model. [169, 170, 183] 

mal expansion (magnetostriction) measurements [169, 170]. The salient feature of the phase 
diagram is a nonmonotonic variation of the Tt(x) curve with a cusp at x^i ~ 0.02 and associated 
with it a splitting in two SC transitions Tci(x) and Tc2(x) in the range x^i < x < Xc2 with Xc2 ~ 
0.045 [184]. This leads to three distinct SC regions A,B and C in the x-T plane. When crossing 
from B to C, Hd exhibits a sudden change in slope [185] and small magnetic moments jU ~ 
10~^/iB seen in /iSi? experiments [167] appear which are absent below Xd (A) and above Xc2- 
There are presently two distinct scenarios, denoted by I and II, to explain these observations. 

(I) The Tc-crossing model 

The most widely accepted interpretation has been given quite early in [5, 186] in terms of a 
Ginzburg-Landau model involving only SC order parameters. The fundamental assumption is 
that two SC phases with gap functions belonging to different representations, e.g. Fi and Fs of 
the cubic symmetry group Oh, have close transition temperatures Tt(Fi,x) and Tt(F5,x) which 
cross at a critical concentration Xd. Below Xd Tc(F5,x) > T(;(Fi,x) and the unconventional F5 
(A) phase which has gap nodes is stable. Above Xd the inequality is reversed and the conven- 
tional (fully symmetric) Fi (B) phase is stable immediately below Td = Tc(Fi,x). For lower 
temperatures (T < Tf2(x)) higher order mixing terms in the GL free energy functional stabilise 
a mixed Fi F5 (C) SC phase with a different Hd, hence the observed kink. The C- phase gap 
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function may be non-unitary and time reversal symmetry breaking, thus leading to condensate 
magnetic moments induced around impurities and defects which allegedly have been observed 
in the juSR experiments. For even larger concentrations enhanced paibreaking or the reduc- 
tion of the spin excitation energy scale Tw(x) involved in pair formation leads to a decrease 
in Tc(ri,x) = Tci and a second crossing point appears. However a GL expansion around x^i 
cannot naturally explain the second crossing. 

(II) The SC/SDW coexistence model 

An alternative scenario to explain the phase diagram in fig. 25 (right panel) has appeared as 
a result of investigations of the differential thermal expansion a = l~^dl/dT [169, 170]. In 
this explanation the SC phases A and B are assumed as identical unconventional SC phases as 
infered from the almost equal pressure dependence away from the critical i . In phase C the SC 
order parameter is unchanged by a SDW modulation which appears within the unconventional 
SC phase. This was previously already proposed theoretically in [171] within the context of 
a model calculation for coexistence of a necessarily unconventional SC and a SDW phase. 
The moments seen in the C phase are then associated with the SDW. The main reason for this 
alternative interpretation of the x-T phase diagram is the discovery of a new line of anomaly 
Ti(x) in a(T) for subcritical x < x^i which starts at a Ti(0) identical to the temperature for 
which B*(T) = as seen in the B-T phase diagram. With growing x the a(T) anomaly becomes 
increasingly sharp and Tz,(x) exactly joins the Td line at Xc2. This suggests to interpret Tl(x) 
as the onset temperature of short range magnetic correlations which turn into the long range 
SDW state at x^. The a(T) anomaly line is seen to continue even far above Xc2- The shape 
of the OC(T) anomaly above x^i is indeed familiar from other magnetic transitions. It is also in 
accord with the interpretation of a huge (Cr-type) anomaly in the ultrasonic attenuation [187]. 
The SC/SDW coexistence scenario lacks however a natural explanation for the observed kink 
in H^-i at T^z- 

Recent flux-creep experiments [166, 188] have further complicated the picture. A steep 
drop in the flux creep rate was observed below Tf2 which may naturally be explained within the 
Tc- crossing model because the nonunitary C-phase may have domain walls which pin vortices 
very effectively [189]. Using this effect for mapping out the C/A-phase boundary it was found 
that contrary to the original view shown in fig. 25 the C-phase boundary is not perpendicular 
at Xci but rather has a low temperature tail extending down to x = 0, i.e. the drop in flux creep 
indicating a transition into the C-phase is even seen in pure UBeo. 

5.2 NODAL SUPERCONDUCTIVITY IN UPd2Al3 MEDIATED BY MAGNETIC EXCI- 
TONS 

Since its discovery in 1991 by Geibel et al. [163] the extensive experimental studies on this 
antiferromagnetically ordered moderate HE superconductor (see table 4) have lead to two im- 
portant conclusions: Firstly the 5f electrons have a dual nature, partly localised and partly 
itinerant where the former carry the magnetic moments and the latter become superconducting 
[190, 191]. Secondly there is now compelling evidence that the internal low energy excitations 
of the localised 5f subsystem mediate an effective pairing potential between the itinerant elec- 
trons. This leads to a superconducting state which has an anisotropic gap function with line 
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nodes [157]. Moreover the virtual low energy excitations of localised 5f electrons lead to a 
quasiparticle mass enhancement factor m*/mi, which is thought to be the origin of the normal 
HF state of UPd2Al3 [29]. Indeed this new mechanism gives consistent and mass enhance- 
ment factor within a strong coupling approach based on a dual 5f electron model Hamiltonian 
[159]. As discussed in sect. 4, HF superconductivity is assumed to exhibit a non-phononic 
SC pairing mechanism leading to nodal gap functions because the strong on-site repulsion of 
quasi-particles prevents the (nearly) isotropic SC state most favored by the electron-phonon 
mechanism. However, UPd2Al3 is the only case known so far where this theoretical conjecture 
has actually been proven by experiment. This was achieved by complementary INS [192-195] 
and quasiparticle tunneling experiments [196] which, respectively, probe the two aspects of the 
dual nature of 5f electrons: (1) The internal singlet-singlet CEF excitations (5^6 meV) of 
localised 5f electrons which form propagating magnetic exciton bands in the range of 1-8 meV 
due to inter-site interactions. (2) The tunneling current which probes the superconducting gap 
of itinerant 5f electrons in epitaxially grown UPd2Al3-A10jc-Pb tunneling devices. There is a 
considerable interaction between localised and itinerant 5f electrons signified on one hand by 
the appearance of a resonance in the INS scattering function associated with the SC gap. On 
the other hand, even more compelling is the presence of typical 'strong coupling' signatures in 
the tunneling DOS at about the magnetic exciton energy in the center of the AF BZ as deter- 
mined in INS. This is a direct proof that the exchange of these magnetic bosons mediates the 
SC pairing in UPd2Al3. This important new mechanism is distinctly different from both the 
electron-phonon and the spin fluctuation mechanism (sect.4). 

5.2.1 PHYSICAL PROPERTIES 

Further experimental evidence for the dual nature of 5f electrons comes from susecptibility 
measurements [197] above Tc where x(^) for field perpendicular to c shows a pronounced 
maximum below T = 50 K, typical for the effect of CEF states which originate from the 5f^ 
configuration of the ions. The presence of localised states is also seen in Knight shift 
measurements [191] and optical experiments [198]. There in addition the formation of heavy 
itinerant quasiparticles below T* was concluded from an analysis of the Drude peak in the 
optical conductivity. The low lying CEF states of localised 5f electrons are of singlet-singlet 
(doublet) type with a splitting 5^6 meV as obtained from a fit to the overall magnetic exciton 
dispersion [158]. This means that the AF ordering (table 4) found below is of the induced 
moment type. The FM ordered ab planes are stacked along c with an AF wave vector Q = 
(0,0,0.5) (r.l.u.) [199, 200] and an easy axis [100]. Moment reorientation can only be observed 
for fields in the ab plane and the phase diagram was determined in [200]. 

The superconducting state below T^ was investigated in numerous experiments, but a defi- 
nite conclusion on the symmetry of the pair state has not been achieved yet. Upper critical field 
measurements show a flattening of Hc2(T) for low temperatures [201] which is interpreted as 
the effect of Pauli limiting in a spin singlet pair state. This hypothesis was further investigated 
by ^^Al Knight shift experiments [202] which shows a considerable reduction below T^ again 
in favor of singlet pairing, however, the interpretation is not unambiguous due to a large local 
moment contributions to K^. The ^^Al NMR relaxation rate Tj^^ was found to exhibit clear T^ 
power law [202], naively interpreted as evidence for line nodes in the gap function A(k) leading 
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mj 

-molK^ - 




Ta,[K] 


Ai [Mb] 


Tc [K] 




Rlnl 


UPd2Al3 




140 


14.3 


0.85 


1.8 


1.2 


0.67 


UNiaAlg 




120 


4.6 


0.20 


1.2 


0.2-0.4 


0.12 



Table 4: Material parameters of UPd2Al3 and its isostructural sister compound UNiiAls. 



to a quasiparticle DOS ~ |£'|. This was also suggested by "•'Pd NMR/NQR [203]. Conclusions 
from the low temperature specific heat C(T) which should be ~ T^ for line nodes are hampered 
by the difficult subtraction of the nuclear contribution. A C(T) = /qT + aT^ behaviour was found 
where the residual value % scales with the width of the SC transition [204] that characterises 
sample quality. The isostructural sister compound UNiiAla is also a superconductor with co- 
existing incommensurate magnetic order (table 4). However, contrary to UPd2Al3 its 5f states 
have all a delocalised character as is obvious from the much smaller moment, the incommmen- 
surate ordering wave vector and the small entropy release AS at T;v which is only one sixth of 
the value in UPd2Al3. Therefore, superconductivity in this compound is thought to be mediated 
by spin fluctuations of conduction electrons rather than localised 5f excitations as in UPd2Al3. 
However, as in UPt3, a spin triplet gap function has been proposed for UNi2Al3 [205] which 
according to NMR results should have a line node. Note that this is in conflict with the simple 
theory of isotropic AF spin fluctuations proposed for U-HF compounds since they would rather 
prefer spin singlet pairing (sect. 4) [145]. Also it suggests that Blount's theorem which predicts 
the absence of line nodes for strong spin-orbit coupling cannot naively be applied to real U-HF 
superconductors (sect. 4). Finally we note that from recent Ht.2 -measurements on epitaxially 
grown thin films [206] a considerable paramagnetic Pauli-limiting effect was deduced for both 
a and c directions which led the authors to the opposite conclusion that UNi2Al3 should also 
have a SC singlet pair state. 



5.2.2 ELECTRONIC STRUCTURE, THE DUAL MODEL 

Within conventional LDA type band structure calculations for UPd2Al3 [207-209] there is no 
way to treat the partly localised and partly itinerant character of 5f electrons properly, all 5f 
orbitals are incorporated in the basis set for the band states. While the FS topology corresponds 
reasonably to dHvA results [210] the effective masses of the various sheets are far too small. 
Even within a modified LSDA (local spin density approximation) treatment including self in- 
teraction corrections [211] the effective masses are still too small by an overall factor of ten. 
The dual model opens an attractive way to remove this discrepancy [29]. The coupling of itin- 
erant 5f- electrons to the low lying discrete CEF excitations of the localised subsystem offers 
an effective mechanism of mass enhancement (sect. 2). In this approach the calculaction of 
heavy bands proceeds in three steps: 1) 5f orbitals in the jj-coupling limit (j = |) are used and 
those with = ±|,±^ are excluded from the basis set of LDA bands which comprises only jj. = 
±1 states. The band center is fixed to obain the correct 5f count. 2) The multiplet structure of 
localised states 5f is calculated using the proper intra-shell Coulomb interactions. The ground 
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Figure 26: Left panel: Fermi surface of UPd2Al3 calculated within the dual model [29]. The 
main cylinder part has a heavy mass with m* = 19 - 33 m and the torus ('crown') has m* = 65 m. 
Right panel: Comparison of experimental dHvA frequencies (black symbols) from [209] and 
calculated frequencies (open symbols) from the dual model [29]. Large parabolas correspond 
to the corrugated main FS cylinder. Small arc on top of F-A corresponds to the crown sheet. 



state is a doublet |7 — A,J^ — ±3) which is further split by the CEF potential into two singlets 
|r3,4) = ^i\J = = 3) ± |/ = 4;/, = -3)) (53) 

The splitting energy 5 ~ 6 meV is taken as an empirical parameter obtained from INS. 3) The 
scattering of localised band states from CEF excitations leads to the mass renormalisation of 
the former. In the simplest case, without including the dispersion which turns the localised 
CEF transitions with energy 5 into a band of magnetic excitons and neglecting strong coupling 
effects the mass enhancement would be given by 

^^ = 1+2^^ (54, 
o 

* 

with I defined in eqs. (47) and (50). A more refined treatment for ^ appopriate for UPd2Al3 is 

described below. In this approach the different treatment of 5f orbitals with j=| and different jj- 
cannot be explained on the single particle level because their hybridisation matrix elements are 
rather similar. The amplification of the orbital dependence of hybridisation due to many body 
effects has been proposed in [30] and is described in detail in sect. 2.2. 
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Although the starting point is indeed very different from standard LDA, with only a single 
delocalised 5f electron, the main FS sheets are reproduced quite well as is obvious from fig. 26. 
The most prominent FS sheet has the form of a slightly corrugated cylinder oriented along c. 
It has also a large mass enhancement of m*/m ~ 33 (m = free electron mass). In the following 
model discussion one may therefore restrict to this main FS sheet. 

5.2.3 MASS ENHANCEMENT AND SC PAIRING DUE TO MAGNETIC EXCITONS 

The great advantage of the dual model as compared to the pure LDA approach is its natural 
explanation of the mass enhancement factor of m*/m^ ~ 10 (m^ = LDA band mass) which is 
due to the contribution of virtual singlet-singlet CEF excitations to the conduction electron self 
energy. Assuming the value of 5 = 6 meV from INS the proper mass enhancement is predicted 
without further adjustable parameters. In addition to the mass enhancement the interaction 
between the subsystems leads to induced AF magnetic order [157, 158] of localised 5f moments 
and superconductivity in the itinerant part [157-159]. In its most rudimentary form a model 
Hamiltonian is given by eq. (50). In UPd2Al3 the conduction electrons corresponding to the 
main FS cylinder in fig. 26 have only dispersion _L c with e±^^ = e^(kx/ko)^. The magnetic 
exciton dispersion ft)£(q) which originates from the Ft, ^ r4 CEF transitions given by eq. (48) 
is of central importance. For large T eq. (48) describes isolated CEF transitions at an energy 
5 which develop a dispersion at lower temperature due to an effective exchange J(q). In this 
RPA expression a complete softening at the Neel temperature Tj^ in eq. (49) for induced AF 
order is expected at the AF ordering wave vector Q = (0,0, j). Note that, unlike spin waves the 
magnetic excitons which originate in singlet- singlet CEF transitions may exist already above T^y 
and their softening signifies the approaching induced AF order. The experimental investigation 
of magnetic excitons was undertaken in many INS studies [157, 194, 212, 213]. In the former 
the overall dispersion up to the maximum ~ 8 meV was determined and theoretically analysed 
in [158]. The latter focused on high resolution analysis of the low energy excitations around 
the AF wave vector Q. The result is shown in the left panel of fig. 27. Indeed, contrary to the 
above simple RPA singlet-singlet result a magnetic exciton gap of about 1 meV appears. This 
may be caused by the effect of higher magnetic CEF states whose contribution to the staggered 
susceptibility leads to AF order slightly before complete softening is achieved and also by self 
energy effects beyond RPA. The dispersion along c is much larger than in the hexagonal ab- 
plane which may lead to a considerable DOS for the low energy magnetic excitons around 1 
meV. The spectral shape of magnetic excitons is Lorentzian above T^ and evolves into a double 
peak structure where the lower has the appeareance of a sharp resonance within the SC gap 
region [157, 213]. This is a signature of the strong residual interaction between itinerant (SC) 
electrons and the excitations within the localised 5f subsystem. A model analysis of the INS 
spectra at Q leads to conclude that A/(Oe{Q) ~ 1, i.e. SC gap amplitude and magnetic exciton 
gap are nearly degenerate. 

Complementary evidence for this interaction was found in the quasiparticle tunneling spec- 
tra of c-axis oriented epitaxially grown UPd2Al3 films [196] which probe the SC electrons. The 
differential conductivity dl/dV of the UPd2Al3-A10x-Pb tunneling barrier is proportional to the 
quasiparticle DOS Ns(tt)). In the weak coupling BCS limit it shows the monotonic square root 
singular behaviour above A(k|| c). However, retardation effects in the strong coupling limit for 
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0.5 ^z(''-l-U-) 0.6 ~' voii^vj 



Figure 27: Left panel: Magnetic exciton dispersion C!)£(q) (q =(0,0,qj) derived from low en- 
ergy INS [157] for temperatures above (crosses) and below (circles) T^. Right panel: Differ- 
ential conductivity from tunneling for various temperatures (the arrow indicates the sequence 
of curves). The hump at 1 meV is a strong coupling signature of the magnetic exciton at Q = 
(0,0,0.5) (cf. left panel) which mediates Cooper pairing [196]. 



the effective pairing interaction lead to characteristic signatures in dl/dV ~ N5(CQ) at an energy 
which corresponds to typical energies of the bosons which mediate the interaction. This is well 
know in the case of strong-coupling electron-phonon superconductors. In their breakthrough 
tunneling experiment in UPd2Al3 Jourdan et al. [196] have for the first time seen evidence 
of the boson that mediates superconductivity in a HF metal. Its signature can clearly be seen 
in the tunneling spectra in the right panel of fig. 27 at an energy slightly above 1 meV. This is 
far too small for characteristic phonon frequencies, since the Debye energy is of the order 13 
meV. Instead it corresponds directly to the magnetic exciton energy at the AF vector Q for T < 
Tc in the left panel. This strongly suggests that magnetic excitons mediate SC pair formation 
in UPd2Al3. The coupling between the itinerant heavy quasiparticles and the magnetic exciton 
in the center of the AF BZ (q^ = 0.5) shows substantial retardation because the group velocity 
(dcOE I dq^ of the latter is much smaller than the Fermi velocity v^ of heavy quasiparticles due 
to the flatness of 0)E{qz) over a considerable part of the BZ along c* [213]. Within Eliashberg 
theory the quasiparticle DOS can be explained by using a phenomenological retarded potential 
centered at (Oe{Q), the analysis [157] leads to the conclusion 2A/Tc = 5.6 in agreement with the 
gap estimate above and also with NMR results [214]. Together the two complementary results 
shown in fig. 27 present the first direct evidence for a non-phononic SC pairing mechanism in a 
HF superconductor. The magnetic exciton-mediated pairing identified here is distinctly differ- 
ent from both the conventional electron-phonon mechanism and the common spin-fluctuation 
mechanism (sect. 4) which does not involve any localised electron component. 
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5.2.4 SYMMETRY OF THE SC GAP FUNCTION AND ELIASHBERG THEORY 



The question of the symmetry of the gap function A(k) is most difficult to resolve in uncon- 
ventional superconductors, and perhaps with the exception of UPts (see [6, 9] for a review) 
has not been unambiguously achieved for any of the HF superconductors. And yet this is very 
important because the nodal structure of the gap function and the associated low temperature 
thermodynamic and transport properties are directly determined by the symmetry class of the 
gap function. It was suggested [195, 215] that the INS double peak structure for T <S T^ men- 
tioned previously requires the translational symmetry property A(k + Q) = -A(k). Furthermore 
Knight shift and Hci results mentioned above have lead to the proposal of a spin singlet gap 
with A(— k) = A(k). Together this suggest a gap function with node lines A(k) = perpendic- 
ular to c at the AF Bragg planes k = ±^Q (Ai^ in table 5). To investigate the consistency of 
the magnetic exciton model strong coupling calculations based on the Hamiltonian in eq. (50) 
have been performed [159]. Such model calculations cannot predict reliably the symmetry of 
the gap function, but they may decide which one is the most favorable within a restricted class 
of possible gap functions . In addition one may at the same time obtain the mass enhancement 
(m*/mt,) of normal quasiparticles from the self energy and Tc from the Eliashberg equations. 
Both quantities are determined by the effective retarded potential due to magnetic exchange 
given in eq. (51). For simplicity we use an empirical form for (Oe(<i) with a finite gap and no 
dispersion ± c, given by 



The parameters co^, a and /3 (fig. 28) are chosen to describe the experimental dispersion along c 
obtained from INS [212]. The effective interaction of eq. (51) breaks spin rotational symmetry 
in a maximal (Ising type) manner. As explained in sect. 4 this leads to a classification into ESP 
and OSP pair states with spin projection factors p = ±1. The Eliashberg equation for the gap 
function is then given by 



where G(k,i(a^,) is the renormalised conduction electron Green's function and T^ is obtained 
from A(Tc) = 1. The result for the T^ of various gap function candidates as function of the 
dimensionless coupling parameter g are shown in fig. 28. Likewise the dependence of m*/m^ 
on g which is roughly linear may be calculated. From dHvA results the enhancement factor for 
the main FS cylinder in fig. 26 is m*/m^ ~ 10. This would imply g - 2. The corresponding 
theoretical value T^ = 2.9 K for the degenerate Aig and Ai„ states from fig. 28 is a factor 1.6 
larger than the experimental value of 1.8 K. Given that the latter is reduced due to the action 
of static AF order, this result shows that the magnetic exciton mechanism gives both mass 
enhancement and a T^ which are consistent. This supports the previous empirical conjecture for 
this new mechanism of HF superconductivity. 

The dispersion of (Osik) and hence the k-dependence of the interaction is strongest along 
c. Therefore one may restrict to SC candidate states of the type A(k) = A ^(k^) with a few 
possible form factors ^'(k^) given in table 5. It is seen in fig. 28 that the A^, and Ai^, OSP spin 
pairing states have the highest and equal T^'s. This degeneracy between even and odd parity 
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A(r)A(k,ia^) =p^ £ V(k-k',ia)„-ia)0|G(k',i(o;;)|2A(k',ia);). 



(56) 
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Figure 28: Left panel: Critical temperature as function of the dimensionless coupling constant g 
= g/e± = ^^^^cl^- Model parameters for cosiq) are O)^ = 0.01e_L (e_L= conduction band width) 

and j8 = 0.8. Full (dashed) curve: OSP (ESP) SC states (see table 5) [159]. Right panel: DOS 
in the Aig state of table 5 with N(E) ~ E for E < Aq. Inset shows the DOS for field of IT 
and orientation || c (0 = 0). Finite N(0) is induced by the Doppler shift of quasiparticles and 
contributes to specific heat and transport [216]. 



states is due to the Ising type anisotropy of the effective interaction in spin space signified by 
their equal spin projection factors p = -1 (table 5). It will be lifted by the interaction of the SC 
order parameters with the background AF order. The respective node lines for 4>(k^) would be 
at = ±7u/c (zone boundary) for Aig and k^ = for Ai„ (zone center). Note that both states 
have S^ = and therefore both should show Pauli limiting of Hc2 and a Knight shift reduction 
below Tp. For the odd parity state this is due to the fact that the other 'triplet' components 
with S^ = ±1 are not stable, i.e. the d - vector is pinned along c. Note that although Ai^ with 
A(k) = Aqcos/c^ has a node line, it is not an unconventional order parameter in the strict sense 
since it has the full Dg/, symmetry. The quasiparticle DOS of the Aig state is shown in fig. 28 
is linear for E <^ Aq due to the node line at k^ = The inset shows that an external field 
induces a finite residual DOS at E = according to the Doppler shift expression of eq. (42). 
The residual DOS depends on the field direction relative to the nodal position. This leads to a 
polar field- angle dependent residual y(0) -coefficient and thermal conductivity Kii(6). There is 
no azimuthal -dependence due to cylindrical FS symmetry. The residual jq, normalised to the 
normal state Jc„ is given by 

K{(oM) = -[cD^-A{^f]^&H{G)^-A{^^)) (57) 
CO 
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Table 5: Spin and orbital structure of the possible gap functions which are solutions of the 
Eliashberg equations for the dual model of UPd2Al3. The Aig OSP state is fully symmetric and 
its node line (fc^ = ±;r/c) is not enforced by symmetry. 



p 




Dg/, repres. 


spin pairing 




-1 


|o,o) = ^(in)-UT)) 




OSP 


cos(c^j.) 


-1 


|i,o) = ;^(in>+UT» 


rr (Ai„) 


OSP 


sin(c^^) 


+1 


|i,±i) = |TT),|U) 


rr (Ai„) 


ESP 


sin(2c^^) 



Here (» = © — Vs(r, 0) • k is the Doppler shifted energy of the quasiparticle, k and r are its 
momentum and its position with respect to the vortex core oriented along Q. The double av- 
erage is performed over the FS and the inter-vortex region and ®h is the Heaviside function. 
Furthermore f (to) is the effective quasiparticle lifetime. As already mentioned in sect. 3.4 for 
CeCoIn5 angle-resolved measurements of thermal conductivity give important information the 
position of the node lines of A(k) with respect to the crystal axes. This is of great help in the 
determination of gap symmetry. One should have H <S Hc2 to minimise the background ef- 
fect of the Hc2 anisotropy. These experiments, which were proposed in [158], have recently 
been performed on UPd2Al3 single crystals [217] in a rotating field geometry with heat cur- 
rent perpendicular to the field (vortex) direction to probe the itinerant quasiparticle contribution 
(fig. 29). For fields rotated in the hexagonal ab-plane no oscillations have been found excluding 
the presence of node lines perpendicular to the ab plane as has been proposed in the model of 
[218]. On the other hand field rotation in the ac-plane clearly leads to twofold oscillations for 
H ^ Ht2 proving the existence of a node line which is parallel to the ab-plane as is predicted 
by the magnetic exciton model of SC in UPd2Al3 [159] (fig. 29). Because of the cylindrical FS 
geometry the oscillations are qualitatively similar for the three gap functions of table 5. Thus 
thermal transport measurements alone are in this case not able to fix the position of the horizon- 
tal node line. However if one accepts that the Knight shift results point to singlet (OSP) pairing, 
then the Ai^-gap function in table 5 is the proper one for UPd2Al3. 



6 SUMMARY AND OUTLOOK 

The understanding of the heavy fermion state in Ce- and U- intermetallics and its SC phases 
has made great progress in recent years. In the normal state it has become clear that the origin 
of quasiparticle mass enhancement may be profoundly different. It is described by the Kondo 
lattice mechanism in 4f-Ce compounds and by a dual (multiorbital) model of localised and itin- 
erant 5f electrons in U compounds. For Ce sytems the heavy electron bands may be described 
by the renormalised band theory where the single 4f electron is treated as itinerant with a res- 
onant phase shift. Within the dual model for U-compounds, two electrons are treated as being 
localised in CEF states. The remaining 5f electron forms itinerant conduction bands and its 
coupling to magnetic excitons, i.e. propagating CEF excitations, leads to the mass enhance- 
ment. These different approaches lead to realistic Fermi surfaces in cases like CeCu2Si2, UPt3, 
UPd2Al3 and the mass enhancement may be incorporated naturally. 
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Figure 29: Left panel: Angular variation of the normalised b-axis thermal conductivity 
JCyy(H, 0). Field is rotated in the ac-plane ± to the heat current q [217]. Right panel: Theoreti- 
cal calculation of K"yj(H, 0) according to eq. (57) for the corrugated FS cylinder of UPd2Al3 and 
the Aig gap function. Twofold oscillations in 6 proves the existence of a node line in the ab 
plane [216]. 

The Kondo lattice mechanism in Ce-compounds implies the competition between HF state 
and AF ordering as illustrated in the Doniach phase diagram. When AF order is destroyed by 
variaton of a control parameter, pressure or chemical substitution, the normal state behaviour 
around the QCP shows distinct NFL behaviour. It is characterised by anomalous scaling expo- 
nents in the temperature and field dependences of specific heat, resistivity, thermal expansion 
and other properties. In Ce-compounds these NFL exponents may largely be understood by 
considering the dressing of quasiparticle with a spectrum of soft magnetic spin fluctuations, 
this picture seems especially appropriate for CeNi2Ge2 but also for CeCu2Si2 and CePd2Si2. 
Simultaneously effective quasiparticle interactions mediated by soft spin fluctuations favor SC 
pair formation as predicted by Eliashberg theory. In fact in many Ce-compounds SC domes 
exist around the QCP and frequently but not necessarily is associated with NFL behaviour. 

Mass enhancement and SC pairing are closely related. The importance of virtual magnetic 
exciton exchange for SC pair formation in HF-U compounds with dual 5f electron behaviour 
may therefore be supposed. Indeed this novel mechanism has been identified in the model 
compound UPd2Al3 by complementary tunneling and INS measurements. The former ex- 
hibit a characteristic strong coupling feature at the energy of the zone center magnetic exciton 
mode which identifies the latter as the exchanged boson that mediates pairing. This obser- 
vation is quite unique for HF-compounds, it has not been achieved for the spin fluctutation 
Ce-compounds which is explicable since the latter are overdamped modes whereas magnetic 
excitons are propagating modes with a real group velocity. Whether SC in other U-HF com- 
pounds may be explained with a similar mechanism is a matter of debate. UBei3 is the only 
SC U compound which possibly exhibits a magnetic QCP and may therefore be more similar 
to Ce-based SC, while UPts and URu2Si2 where SC is embedded in a magnetic (or hidden 
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order) phase are more plausible candidates for the magnetic exciton mechanism. Also the AF 
spin fluctuation model in its simple versions strongly favors singlet pairing whereas a number 
of U-HF compounds, notably UPts, exhibit triplet pairing. The latter is not disfavored in the 
magnetic exciton model due to an inherent spin space anisotropy of the effective interaction. 

The most important property of unconventional HF SC is the symmetry class of the gap 
function. Progress in understanding these symmetries for specific compounds has been slow in 
the past since the traditional method, i.e. interpreting 'power law' behaviour of thermodynamic 
and transport quantities etc. is ambiguous and often contradictory. Recently the situation has 
improved considerably with the advent of angle-resolved specific heat and magnetotransport 
measurements that can identify or restrict nodal positions of the gap functions in k-space. To- 
gether with Knight shift measurements they are the most powerful tool to investigate the gap 
symmetry in ultra-pure stoichiometric samples. Up to now the SC gap is known with some 
certainty only for three cases: f-wave triplet (E2„) in UPts, nodal singlet (Aig) in UPd2Al3 and 
d-wave singlet (Bi^ or ) in CeCoIns. 

An issue of permanent interest in HF SC is the coexistence/competition with long range AF 
or SDW order away from the QCP. This is now fairly well understood for the 'A-phase' SDW 
state in CeCu2Si2 and CeCu2(Sii_;,:Ge;,^)2, both experimentally and within theoretical models 
involving a realistic FS geometry. A similar approach may be possible for Cel 15 materials like 
CeRhi_;cCoj(;In5. 

The vortex phase can give further insight into the SC state. For UPt3 the appearance of three 
distinct SC regions in the B-T plane was a sure sign of a multicomponent SC order parameter. 
More recently the observation of a distinct SC vortex phase in the low temperature and high field 
comer close to Hc2 for CeCoIns has been interpreted as the long-sought FFLO phase where SC 
pairs have finite momentum wich leads to a segmentation of vortices. 

Many fundamental issues in HF materials, like the importance of unconventional density 
wave states for NFL behaviour, QCP scenarios with a local breakup of quasiparticles, the effect 
of time reversal and inversion symmetry breaking in unconventional SC states, the vortex and 
FFLO phase for nodal SC gap functions and many others, are still incompletely understood. 
This guarantees that heavy fermion physics will be a thriving field in the future. 
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List of acronyms 



AF 


antiferromagnet, antiferromagnetism, antiferromagnetic 


AFBZ 


antiferromagnetic Brillouin zone 


ARPES 


angle resolved photoemission spectroscopy 


BIS 


bremsstrahlung isochromat spectroscopy 


BZ 


Brillouin zone 


CDW 


charge density wave 


CEF 


crystalline electric field 


CPT 


cluster perturbation theory 


dHvA 


de Haas- van Alphen 


DLRO 


diagonal long range order 


DOS 


density of states 


ESP 


equal spin pairing 


FFLO 


Fulde-Ferrell-Larkin-Ovchinnikov 


GL 


Ginzburg-Landau 


LFL 


Landau-Fermi liquid 


FLEX 


fluctuation exchange 


EM 


ferromagnet, ferromagnetism, ferromagnetic 


FS 


Fermi surface 


HE 


heavy fermion 


IC 


incommensurate 


INS 


inelastic neutron scattering 


LDA 


local density approximation 


LSDA 


local spin density approximation 


mi 


mean field 


NFL 


non-Fermi liquid 


n.n. 


next neighbor 


n.n.n. 


next nearest neighbor 


NMR 


nuclear magnetic resonance 


NQR 


nuclear quadrupole resonance 


ODLRO 


off-diagonal long range order 


OSP 


opposite spin pairing 


PES 


photoemission spectroscopy 


QCP 


quantum critical point 


RPA 


random phase approximation 


RKKY 


Ruderman-Kittel-Kasuya-Yoshida 


SC 


superconductor, superconductivity, superconducting 


SIC-LDA 


self interaction corrected local density approximation 


SDW 


spin density wave 


s.o. 


spin orbit 


TB 


tight binding 
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